Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1400 1777 (N38W61 L=152)  1483    180        0 (1326/1719-1479/1810)  -4.8 W49-W73 N34-N39
      2 1351 1654 (N33W45 L=136)   841    134        0 (1302/1594-1394/1713)  71.6 W42-W48 N28-N37
      3 1489 1522 (N21W46 L=137)   265     59        0 (1461/1502-1515/1540) -56.3 W44-W47 N19-N22
      4 1737  992 (S15W44 L=135)   733    193        1 (1675/ 950-1804/1039) -42.0 W39-W50 S19-S10
      5 1860  951 (S22W57 L=148)   761     97        1 (1815/ 926-1886/ 990)   5.7 W51-W61 S24-S20
      6 1575  866 (S17W29 L=120)   740    200        0 (1551/ 772-1589/ 961)  61.4 W25-W31 S22-S11
      7  233  772 (N09E63 L= 28)  1783    266        6 ( 179/ 669- 280/ 884) -70.3 E67-E59 N02-N16
      8  276  819 (N12E56 L= 35)   481     57        2 ( 249/ 793- 305/ 849) -36.9 E59-E53 N10-N13
      9 1408  682 (S23W14 L=105)   275     88        0 (1383/ 640-1423/ 726)  45.0 W11-W15 S25-S20
     10  283  615 (S01E64 L= 27)   487     80        2 ( 255/ 589- 312/ 654)  11.3 E68-E58 S02-N01
     11  908  404 (S27E26 L= 65)   521    112        1 ( 863/ 372- 947/ 439)  14.0 E29-E21 S28-S25
     12 1328  289 (S50E03 L= 88)   952    250        0 (1229/ 225-1424/ 353)  10.8 E14-W07 S53-S47
  

    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