Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1190 1495 (N25W06 L=319)   361    134        0 (1166/1427-1205/1552) -69.4 W03-W07 N20-N29
      2  447 1393 (N12E42 L=272)  1026    239        4 ( 398/1328- 503/1478) -42.0 E47-E36 N08-N17
      3  470 1278 (N05E38 L=276)   296     81        1 ( 442/1247- 504/1318)  53.1 E39-E35 N03-N07
      4 1659 1186 (N12W39 L=353)  1803    436        4 (1544/1079-1806/1266) -15.9 W30-W52 N06-N16
      5  611  984 (S11E25 L=288)   255     85        0 ( 567/ 965- 643/1007)  36.9 E27-E22 S12-S09
      6  488  905 (S17E33 L=280)   215     47        0 ( 459/ 889- 511/ 920) -14.0 E35-E31 S17-S15
      7 1566  826 (S11W37 L=350)   423    105        0 (1524/ 802-1620/ 852)   0.0 W33-W40 S11-S09
      8 1117  643 (S27W10 L=324)  2160    620        4 ( 820/ 531-1386/ 809)  34.9 E09-W24 S37-S14
      9  350  632 (S35E47 L=267)   337     54        0 ( 327/ 607- 373/ 654) -45.0 E49-E45 S36-S33
  

    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