Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  574 1423 (N12E37 L=103)   365     94        0 ( 548/1393- 610/1461) -26.6 E40-E34 N10-N14
      2 1866 1257 (N32W59 L=199)   655    118        2 (1838/1237-1894/1282)  -9.5 W55-W62 N30-N32
      3  462 1222 (S02E40 L=101)   570    151        1 ( 422/1144- 511/1278) -29.1 E43-E34 S05-N01
      4 1711 1155 (N21W43 L=183)  1084    110        3 (1665/1083-1750/1217) -45.0 W38-W46 N16-N23
      5  139 1141 (S12E71 L= 70)  1974    167        2 ( 118/1047- 163/1199) -83.7 E76-E66 S17-S08
      6 1171 1163 (N09W06 L=146)   324     98        0 (1107/1138-1213/1185)   0.0 W01-W08 N08-N10
      7  438  938 (S19E36 L=104)   687    223        1 ( 359/ 872- 507/1011)  50.7 E41-E31 S24-S13
      8  758  853 (S17E12 L=128)   222     74        1 ( 732/ 827- 798/ 880)  36.9 E13-E09 S19-S15
      9  480  819 (S25E31 L=109)  1089    381        0 ( 373/ 722- 624/ 869)  35.8 E39-E21 S33-S20
     10 1299  765 (S11W22 L=162)   313     97        0 (1242/ 748-1347/ 788)  33.7 W18-W24 S13-S08
     11 1552  439 (S25W53 L=193)  1714    261        4 (1472/ 358-1631/ 553) -85.9 W49-W56 S30-S16
     12 1092  434 (S37W20 L=160)   733    162        0 (1006/ 420-1172/ 453)  20.0 W14-W25 S39-S33
  

    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