Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  710 1649 (N27E38 L=245)  1070    349        2 ( 618/1559- 775/1746) -23.2 E51-E30 N22-N32
      2  985 1658 (N34E21 L=263)   299     52        1 ( 955/1640-1018/1674)  14.0 E22-E18 N32-N35
      3  484 1587 (N18E52 L=231)   373     59        0 ( 464/1557- 513/1618)  90.0 E53-E50 N15-N20
      4  366 1472 (N08E57 L=227)  1024    143        2 ( 331/1416- 415/1524) -50.2 E60-E51 N05-N11
      5  332 1130 (S13E48 L=236)   354     44        2 ( 311/1099- 354/1157)  90.0 E49-E45 S14-S10
      6  268 1069 (S17E52 L=231)   939    127        5 ( 234/1017- 303/1123)  80.5 E54-E48 S20-S14
      7  200 1040 (S20E59 L=225)   895    218        2 ( 158/ 937- 226/1110)  76.0 E61-E56 S27-S15
      8  819  919 (S14E09 L=274)   358     74        1 ( 783/ 886- 854/ 951)  68.2 E11-E07 S16-S11
      9  275  815 (S33E48 L=235)   352     69        2 ( 250/ 795- 303/ 831)  11.3 E50-E45 S34-S31
     10  773  692 (S29E07 L=276)   323     94        0 ( 729/ 659- 811/ 721)  51.3 E09-E05 S31-S25
     11  608  565 (S42E17 L=266)   563    127        1 ( 557/ 526- 652/ 598)  36.9 E21-E13 S45-S38
  

    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