Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1282 1837 (N67W14 L= 24)   664     76        1 (1247/1809-1316/1855)   0.0 W08-W19 N63-N69
      2 1434 1446 (N33W24 L= 34)  1255    322        2 (1381/1326-1542/1541) -49.8 W19-W32 N25-N39
      3  508 1437 (N20E42 L=328)   400     79        0 ( 471/1414- 549/1460) -15.9 E45-E38 N18-N20
      4  865 1412 (N23E16 L=354)   281     62        1 ( 848/1388- 887/1444)  90.0 E17-E13 N21-N24
      5 1730 1346 (N30W51 L= 61)   406     47        2 (1708/1315-1756/1377)  63.4 W48-W54 N27-N31
      6 1333 1285 (N21W17 L= 27)   321    129        1 (1297/1246-1385/1334) -45.0 W13-W20 N18-N23
      7  675 1272 (N11E25 L=345)   240     51        1 ( 645/1254- 704/1293)   0.0 E27-E23 N10-N12
      8 1363 1197 (N15W19 L= 29)   233     75        1 (1329/1176-1393/1217) -11.3 W16-W21 N14-N16
      9 1501 1035 (N07W30 L= 40)   292     68        1 (1485/1007-1515/1071) -76.0 W29-W31 N05-N09
     10  872  807 (S14E07 L=  3)   442    173        0 ( 817/ 745- 922/ 882)  60.9 E09-E04 S18-S09
     11 1490  825 (S06W32 L= 42)   293     83        1 (1468/ 783-1521/ 852)  59.0 W31-W34 S08-S03
     12  462  732 (S25E36 L=334)   944    171        1 ( 437/ 659- 486/ 822)  74.7 E38-E33 S29-S18
     13 1494  711 (S13W35 L= 45)   251     91        0 (1456/ 683-1535/ 736)  45.0 W32-W37 S14-S10
     14 1187  690 (S18W15 L= 25)   292     79        0 (1161/ 657-1223/ 717)  53.1 W13-W16 S19-S15
     15 1426  561 (S23W34 L= 44)   283    120        0 (1411/ 528-1459/ 603) -78.7 W32-W37 S25-S20
     16 1080  550 (S28W10 L= 20)   854    194        1 ( 983/ 521-1156/ 593)  28.6 W03-W14 S31-S24
     17  540  536 (S37E31 L=339)   968    148        2 ( 512/ 472- 584/ 577) -90.0 E33-E26 S42-S34
     18 1536  486 (S27W47 L= 57)  3016    559        2 (1314/ 433-1779/ 568)  18.4 W30-W72 S33-S18
     19 1277  536 (S27W24 L= 34)   374    115        0 (1231/ 500-1323/ 558)  45.0 W22-W27 S29-S24
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y