Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  550 1798 (N29E66 L=235)   724     93        3 ( 509/1761- 577/1836) -49.4 E74-E61 N25-N32
      2  849 1766 (N33E37 L=264)   571     93        2 ( 822/1729- 891/1805) -74.1 E38-E34 N30-N37
      3 1209 1604 (N32E06 L=295)  1080    277        1 (1084/1558-1298/1668)  47.5 E13-E01 N28-N40
      4 1167 1382 (N17E02 L=300)   278     87        2 (1129/1349-1201/1414)  59.0 E03-W00 N14-N19
      5 1546 1259 (N21W25 L=326)  1041    296        0 (1461/1191-1643/1301)   0.0 W19-W32 N14-N23
      6 1874  933 (N14W61 L=  2)  1288    129        2 (1858/ 862-1894/1001) -82.9 W58-W63 N09-N18
      7 1048  928 (S12W04 L=305)   335    139        0 ( 994/ 893-1113/ 956)   7.1 W00-W08 S12-S10
      8  994  739 (S24W06 L=308)   851    313        2 ( 890/ 681-1097/ 830)  58.4 W01-W10 S29-S15
      9 1193  695 (S20W19 L=321)   339    124        0 (1138/ 659-1234/ 724)   0.0 W15-W22 S21-S19
     10 1556  641 (S12W44 L=345)   571    156        1 (1493/ 602-1631/ 676)  41.2 W40-W49 S16-S09
  

    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