Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  988 1615 (N33E05 L=293)   631     99        1 ( 958/1561-1022/1662)  90.0 E06-E02 N29-N36
      2  575 1625 (N33E36 L=262)   362     66        0 ( 543/1609- 614/1641)   9.5 E39-E32 N31-N34
      3 1417 1483 (N26W25 L=323)   783    225        0 (1332/1443-1517/1541) -22.6 W19-W31 N23-N29
      4 1594 1310 (N16W36 L=334)   901    284        0 (1469/1234-1706/1368) -22.4 W27-W44 N11-N18
      5  225 1326 (N13E60 L=238)   559     95        0 ( 180/1299- 268/1356)  11.3 E65-E55 N11-N14
      6 1888 1207 (N12W64 L=  2)  1914    213        0 (1860/1125-1932/1318) -65.6 W60-W69 N07-N18
      7 1755 1251 (N13W49 L=347)   355     62        0 (1736/1212-1778/1283) -63.4 W47-W51 N11-N15
      8  174 1247 (N08E64 L=234)   615     56        1 ( 135/1224- 210/1267)  -5.2 E70-E59 N06-N09
      9 1587 1190 (N08W35 L=333)   253     97        0 (1536/1173-1630/1217)  23.2 W31-W38 N07-N10
     10 1610 1038 (S01W37 L=335)   599    143        1 (1584/ 965-1643/1111) -60.3 W34-W39 S04-N03
     11 1849  820 (S11W61 L=359)  1597    149        6 (1798/ 750-1905/ 873)  28.6 W56-W67 S15-S07
     12 1004  674 (S25W00 L=298)  1208    332        0 ( 928/ 559-1144/ 750)  47.5 E04-W08 S32-S20
     13 1662  588 (S27W50 L=348)  1063    141        1 (1583/ 561-1724/ 619)   8.7 W42-W55 S28-S24
     14  265  573 (S33E63 L=235)   720     97        0 ( 229/ 540- 309/ 595)   4.1 E71-E57 S34-S31
     15 1039  535 (S34W03 L=301)   367     82        0 ( 992/ 516-1086/ 557)  -9.5 W00-W06 S35-S32
     16 1481  372 (S43W44 L=342)  2784    377        2 (1305/ 347-1662/ 414)   9.5 W26-W65 S45-S38
     17  669  339 (S50E31 L=268)  1913    307        2 ( 504/ 310- 802/ 387)   3.3 E51-E16 S52-S46
  

    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