Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  865 1817 (N52E34 L=223)   475     71        1 ( 826/1799- 904/1835)  38.7 E36-E31 N50-N54
      2  585 1267 (N09E31 L=225)  2077    623        1 ( 441/1054- 727/1495) -35.4 E49-E18 S02-N20
      3  274 1240 (N03E57 L=200)   591    132        1 ( 228/1196- 324/1294)  36.9 E60-E52 N00-N06
      4 1365 1231 (N17W19 L=275)   283    110        0 (1308/1187-1402/1261) -23.2 W14-W21 N15-N19
      5 1093 1212 (N12W02 L=258)   553    195        2 ( 987/1184-1181/1246)   0.0 E05-W07 N11-N13
      6 1239  845 (S08W16 L=272)  2342    641        6 (1130/ 594-1352/1095) -76.9 W06-W24 S23-N07
      7  292 1007 (S11E50 L=206)   858    161        1 ( 274/ 922- 307/1081) -73.3 E52-E48 S16-S06
      8  642  914 (S12E22 L=234)   436     88        1 ( 593/ 896- 693/ 931)  15.9 E25-E18 S12-S10
      9  218  886 (S20E57 L=199)   406    104        0 ( 191/ 850- 250/ 925) -29.7 E60-E53 S21-S17
     10  931  844 (S12E03 L=254)   557    215        1 ( 888/ 755-1019/ 905)  60.9 E05-W01 S17-S07
     11 1096  527 (S30W13 L=270)   562    130        1 (1039/ 502-1169/ 554)  32.0 W08-W17 S32-S27
  

    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