Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1622 1581 (N21W58 L=182)   617    148        1 (1597/1518-1636/1627)  39.8 W54-W61 N18-N23
      2 1796 1143 (S07W54 L=178)  1070    294        1 (1714/1008-1834/1235) -75.1 W49-W57 S15-N00
      3  158 1211 (N35E62 L= 62)   605    176        0 ( 124/1188- 185/1228)  -9.5 E65-E59 N33-N35
      4 1689  697 (S31W36 L=160)  2048    361        3 (1549/ 566-1822/ 842)  19.1 W23-W50 S35-S25
      5 1141  625 (S19E03 L=121)   267    103        0 (1085/ 611-1194/ 638) -18.4 E06-W00 S20-S17
      6  859  504 (S20E23 L=101)  1130    335        3 ( 699/ 448-1007/ 561) -36.9 E31-E15 S26-S13
      7  700  430 (S21E36 L= 87)   868    201        2 ( 635/ 347- 743/ 506) -90.0 E38-E33 S27-S15
      8  442  433 (S16E57 L= 67)   506     99        1 ( 413/ 390- 469/ 471) -90.0 E58-E54 S19-S13
      9  530  400 (S20E51 L= 73)  1524    324        4 ( 472/ 340- 592/ 457)   8.7 E58-E45 S24-S16
  

    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