Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  839 1423 (N18E20 L=101)   887    227        2 ( 780/1343- 899/1504) -41.6 E25-E14 N14-N23
      2  311 1339 (N03E57 L= 65)   567    160        1 ( 260/1292- 386/1370)  31.0 E60-E50 N00-N06
      3  503 1159 (S04E35 L= 86)   560    110        1 ( 448/1140- 549/1178)  33.7 E39-E31 S05-S01
      4 1843 1052 (N17W59 L=181)  1820    372        2 (1752/ 929-1895/1173) -31.4 W48-W68 N10-N23
      5  930  936 (S09E04 L=118)   589    166        2 ( 868/ 880- 987/ 994)  58.0 E06-E01 S12-S04
      6 1232  805 (S10W17 L=138)   519    203        0 (1151/ 750-1288/ 867)  60.9 W13-W18 S14-S05
      7 1372  832 (S06W25 L=146)   238     32        0 (1353/ 810-1397/ 850)  71.6 W23-W26 S07-S04
      8 1595  770 (S05W42 L=163)   274     57        0 (1581/ 743-1613/ 802) -53.1 W40-W43 S06-S02
      9 1257  744 (S13W20 L=141)   298    101        0 (1238/ 697-1284/ 789)  90.0 W19-W20 S16-S10
     10 1197  683 (S18W18 L=139)   280    128        0 (1147/ 654-1235/ 713)  -9.5 W14-W20 S19-S16
     11 1518  615 (S15W42 L=163)   259     88        1 (1479/ 581-1545/ 650)  78.7 W40-W43 S18-S13
  

    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