Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  484 1421 (N37E37 L=130)   321     61        1 ( 451/1401- 524/1434)   0.0 E40-E33 N35-N38
      2  850 1364 (N29E07 L=159)   238     55        0 ( 823/1343- 879/1385)   0.0 E09-E05 N28-N30
      3 1110 1323 (N23W10 L=176)   205     99        1 (1060/1313-1163/1334) -15.9 W06-W13 N21-N24
      4 1649 1207 (N08W44 L=211)   325     81        0 (1633/1167-1666/1247)  45.0 W42-W46 N05-N09
      5  466 1113 (N17E36 L=130)  1133    408        1 ( 407/1000- 543/1194) -71.6 E41-E30 N10-N22
      6  294  906 (N06E53 L=114)   956    135        3 ( 255/ 841- 335/ 955)  37.9 E57-E48 N01-N08
      7 1585  393 (S47W42 L=208)  1382    151        2 (1519/ 358-1650/ 436)   7.6 W33-W49 S50-S44
      8  499  353 (S33E56 L=110)   573    106        1 ( 483/ 305- 514/ 404)  32.7 E65-E51 S37-S28
      9  605  287 (S39E51 L=116)  1246    122        2 ( 551/ 255- 657/ 319)  10.8 E63-E42 S42-S36
  

    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