Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  829 1786 (N62W18 L=244)   741    102        1 ( 803/1738- 864/1825)  19.7 W12-W26 N58-N64
      2 1290 1556 (N29W34 L=260)   359     75        1 (1259/1531-1334/1576)   0.0 W31-W38 N28-N30
      3 1917  977 (S24W67 L=292)   814     66        1 (1900/ 935-1942/1004) -78.7 W63-W72 S26-S21
      4  203  852 (N15E64 L=162)   483     94        0 ( 180/ 806- 221/ 896) -71.6 E66-E62 N11-N17
      5 1859  837 (S30W56 L=282)  1156    141        3 (1817/ 787-1911/ 878)   4.1 W50-W66 S31-S27
      6 1291  768 (S15W09 L=234)   497    123        3 (1252/ 725-1327/ 804)  68.2 W06-W11 S17-S12
      7 1275  575 (S26W03 L=228)   436    173        0 (1212/ 539-1310/ 639)  38.7 E01-W05 S27-S21
      8  856  551 (S16E23 L=203)   389    130        0 ( 804/ 529- 921/ 583) -50.2 E25-E18 S18-S12
      9  804  436 (S22E31 L=195)  1027    215        5 ( 718/ 408- 893/ 478) -14.0 E36-E24 S25-S17
  

    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