Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1495 1720 (N56W41 L=163)  1041    125        2 (1448/1673-1551/1748)  19.1 W30-W56 N51-N60
      2 1023 1625 (N40E09 L=113)   463    100        2 ( 977/1594-1061/1644)   0.0 E12-E05 N37-N42
      3  797 1414 (N22E20 L=102)   541    155        0 ( 771/1349- 815/1492) -60.9 E22-E17 N18-N27
      4 1335 1452 (N32W17 L=139)   239     81        1 (1301/1419-1354/1488) -53.1 W14-W18 N29-N33
      5  577 1403 (N19E35 L= 87)   464     93        1 ( 555/1366- 609/1442) -78.7 E37-E32 N16-N21
      6  701 1357 (N18E25 L= 97)   492    152        1 ( 670/1278- 726/1422)  90.0 E26-E23 N12-N21
      7  594  780 (S19E25 L= 97)   770    194        2 ( 540/ 694- 654/ 831)  48.4 E29-E21 S25-S15
      8 1706  648 (S15W55 L=177)   443     24        1 (1682/ 626-1729/ 671)  -9.5 W52-W58 S16-S13
      9 1287  417 (S35W29 L=151)   383     88        0 (1245/ 394-1327/ 446)  51.3 W26-W31 S37-S32
     10  817  225 (S61E08 L=115)   483     72        0 ( 788/ 192- 835/ 257) -60.3 E10-E05 S65-S57
  

    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