Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  475 1629 (N24E57 L=194)   510     81        0 ( 457/1585- 501/1666) -59.0 E59-E53 N21-N26
      2  395 1553 (N18E59 L=191)  1661    244        2 ( 332/1454- 443/1641) -18.4 E80-E49 N12-N22
      3  318 1443 (N10E61 L=190)   547    110        1 ( 284/1396- 346/1478) -15.3 E66-E55 N07-N11
      4  256 1350 (N03E63 L=187)   584    145        0 ( 197/1317- 317/1383)  22.6 E70-E57 N00-N05
      5 1735 1296 (N30W46 L=297)   956    103        3 (1695/1248-1763/1340) -33.7 W42-W49 N27-N33
      6  758 1279 (N08E21 L=229)   261    102        0 ( 744/1230- 772/1323) -71.6 E22-E20 N04-N10
      7 1823  836 (N04W61 L=311)   511     46        0 (1797/ 816-1851/ 858)   0.0 W57-W65 N03-N05
      8  388  730 (S32E40 L=210)   595    129        1 ( 356/ 698- 451/ 774) -12.5 E43-E33 S34-S29
      9 1178  619 (S22W19 L=269)   346     69        0 (1144/ 601-1216/ 643)  -9.5 W15-W21 S23-S21
     10  817  536 (S36E03 L=247)   249     58        1 ( 802/ 511- 834/ 567) -45.0 E04-E01 S37-S34
  

    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