Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1403 1741 (N36W57 L=157)  2478    352        3 (1267/1677-1557/1778) -11.6 W40-W74 N28-N39
      2 1267 1614 (N33W36 L=136)   783    139        2 (1229/1546-1317/1674)  78.7 W33-W39 N27-N37
      3  641 1149 (N23E20 L= 80)   274     59        0 ( 612/1115- 666/1178) -76.0 E21-E18 N21-N25
      4  210 1018 (N24E60 L= 40)   659    131        2 ( 160/ 973- 267/1046) -24.8 E67-E53 N20-N26
      5 1693  923 (S18W39 L=139)  1091    354        4 (1544/ 872-1815/1001) -35.5 W29-W50 S23-S08
      6 1830  779 (S32W52 L=152)  2891    349        3 (1755/ 657-1939/ 913)  21.0 W42-W72 S38-S26
      7  256  888 (N16E56 L= 43)   316     47        2 ( 235/ 865- 277/ 911)  26.6 E58-E53 N14-N17
      8  211  803 (N11E64 L= 36)  1468    149        3 ( 182/ 744- 242/ 875)  82.9 E69-E59 N07-N15
      9  150  819 (N13E73 L= 26)  1753    150        2 ( 119/ 744- 170/ 874) -69.4 E78-E69 N08-N16
     10 1461  774 (S19W19 L=119)   320    108        1 (1446/ 718-1479/ 821)  51.3 W17-W21 S22-S16
     11  394  578 (S05E54 L= 46)   606    118        0 ( 343/ 561- 445/ 598) -20.6 E58-E48 S06-S02
     12  798  367 (S26E34 L= 65)   505     72        2 ( 760/ 335- 832/ 401)   8.1 E38-E31 S28-S24
  

    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