Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  381 1641 (N34E68 L=306)   672     69        1 ( 369/1608- 401/1677) -22.6 E76-E63 N31-N36
      2  597 1590 (N33E40 L=335)  1898    385        5 ( 463/1449- 703/1657) -20.1 E56-E26 N24-N37
      3  491 1469 (N24E43 L=331)   789    276        1 ( 414/1384- 536/1583) -85.9 E48-E41 N17-N31
      4  413 1481 (N23E51 L=324)   366     32        0 ( 396/1451- 428/1509) -53.1 E53-E49 N21-N25
      5 1050 1416 (N26E02 L= 13)   342    111        1 (1019/1371-1093/1445)  71.6 E03-W00 N22-N28
      6  885 1344 (N20E12 L=  3)   410    193        0 ( 781/1314- 960/1363)  -4.8 E18-E06 N18-N20
      7  929 1060 (N03E06 L=  9)   250     51        1 ( 898/1037- 961/1077)  18.4 E07-E04 N01-N03
      8 1246 1017 (N03W14 L= 28)   630    189        0 (1180/ 965-1330/1051)  18.4 W09-W18 N00-N05
      9  777 1029 (S01E15 L=360)   206     90        0 ( 737/1016- 821/1048)   0.0 E17-E12 S01-N00
     10  153  918 (S14E68 L=307)  1233    149        2 ( 127/ 850- 174/ 979) -54.5 E72-E64 S18-S09
     11 1609  716 (S12W43 L= 57)  1928    419        3 (1483/ 616-1725/ 898)  55.6 W35-W49 S19-N00
     12  317  636 (S30E52 L=323)   908    151        1 ( 273/ 575- 375/ 676)  26.6 E60-E44 S35-S26
  

    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