Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1289 1801 (N58E06 L=217)   808     73        2 (1271/1753-1311/1843) -61.4 E10-E03 N52-N63
      2 1024 1632 (N34E18 L=206)   674    227        1 ( 978/1537-1092/1711) -31.0 E24-E09 N28-N38
      3 1543 1548 (N44W24 L=248)   657    105        1 (1508/1503-1576/1583) -56.3 W20-W27 N40-N46
      4 1552 1326 (N28W27 L=250)   281    122        0 (1536/1270-1580/1379) -50.2 W24-W29 N25-N31
      5 1809 1227 (N30W53 L=276)  1299    228        0 (1735/1137-1880/1301) -18.4 W43-W62 N26-N32
      6  220  770 (S37E60 L=164)   746     75        1 ( 195/ 735- 249/ 802)  32.0 E64-E54 S39-S34
      7  279  579 (S50E59 L=165)   589     58        1 ( 256/ 551- 308/ 604)  22.6 E64-E51 S52-S47
      8 1348  443 (S29W39 L=263)  2435    466        5 (1129/ 380-1542/ 507)  36.9 W26-W51 S37-S19
  

    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