Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1124 1722 (N52W09 L=210)   309     42        1 (1104/1694-1144/1740) -76.0 W06-W11 N49-N53
      2  556 1515 (N34E38 L=164)  1805    419        2 ( 405/1446- 732/1594) -11.9 E59-E21 N29-N38
      3  743  687 (S18E18 L=184)   496    213        0 ( 648/ 664- 838/ 718)  -4.4 E24-E11 S19-S16
      4  401  581 (S27E48 L=154)  2472    317        3 ( 319/ 426- 474/ 689)  29.6 E76-E38 S40-S19
      5 1709  647 (S21W53 L=254)   417     19        0 (1688/ 627-1730/ 669) -31.0 W50-W55 S22-S19
      6  799  608 (S24E15 L=187)   453    171        0 ( 753/ 547- 874/ 651)  37.9 E18-E09 S27-S20
      7 1186  499 (S31W12 L=214)   265     88        0 (1145/ 481-1232/ 518)  18.4 W09-W15 S32-S29
      8 1574  415 (S38W50 L=251)   780    113        0 (1546/ 368-1597/ 468) -63.4 W46-W52 S41-S33
      9 1460  362 (S42W40 L=242)  1165    117        2 (1407/ 338-1524/ 395)   4.4 W34-W47 S44-S39
     10 1226  308 (S47W20 L=221)  1178    211        1 (1139/ 286-1319/ 346)   0.0 W11-W28 S48-S43
  

    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