Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1469 1733 (N33W62 L=137)  1153     98        1 (1419/1697-1515/1765) -71.6 W58-W65 N29-N35
      2  464 1461 (N48E31 L= 44)   304     59        1 ( 445/1432- 482/1485) -76.0 E33-E29 N45-N49
      3 1835 1129 (S11W57 L=133)   582     82        1 (1802/1095-1866/1166) -59.0 W54-W60 S13-S08
      4 1755  981 (S17W45 L=121)  1348    220        2 (1731/ 878-1775/1096)  69.0 W41-W48 S22-S09
      5  606 1025 (N17E25 L= 51)   231     38        2 ( 584/1004- 631/1047) -71.6 E25-E23 N15-N18
      6  349  861 (N13E48 L= 28)   457     61        1 ( 327/ 816- 366/ 895)  78.7 E49-E45 N10-N15
      7  173  745 (N08E73 L=  3)  1643    147        2 ( 125/ 694- 221/ 814) -41.6 E80-E65 N04-N12
      8 1427  309 (S53W07 L= 82)   399     59        0 (1398/ 294-1463/ 323)  -8.1 W03-W10 S54-S51
      9  902  274 (S36E34 L= 42)   593    111        0 ( 855/ 240- 946/ 320)  18.4 E39-E27 S37-S33
  

    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