Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  741 1694 (N38E27 L=355)  2179    326        4 ( 589/1638- 881/1738) -12.1 E43-E15 N33-N40
      2 1303 1647 (N38W16 L= 39)   305     82        0 (1263/1635-1344/1665) -18.4 W13-W19 N37-N39
      3 1287 1536 (N29W14 L= 37)   277     99        0 (1243/1508-1328/1570) -26.6 W11-W17 N27-N31
      4 1095 1366 (N17W02 L= 25)   995    218        1 (1036/1271-1145/1451)  71.6 E01-W05 N10-N22
      5  284 1336 (N11E55 L=328)   380     53        0 ( 263/1309- 307/1365) -33.7 E58-E51 N09-N13
      6  411 1303 (N10E43 L=340)   335     78        2 ( 386/1265- 437/1340) -63.4 E45-E40 N07-N11
      7 1302 1288 (N13W15 L= 38)   637    241        0 (1199/1236-1407/1338) -17.1 W08-W22 N10-N15
      8  510 1282 (N08E35 L=348)   245     97        0 ( 471/1261- 563/1294)   8.1 E37-E30 N07-N09
      9  300 1164 (N01E50 L=333)   375     65        0 ( 273/1133- 331/1188)  26.6 E52-E47 N00-N02
     10  566 1109 (S02E29 L=354)   484    126        0 ( 516/1052- 614/1154)  45.0 E32-E25 S05-N01
     11  286 1013 (S08E51 L=332)   547     80        1 ( 271/ 962- 311/1053)  80.5 E52-E48 S11-S05
     12  203  960 (S11E60 L=323)   440     41        0 ( 186/ 931- 219/ 994)  56.3 E61-E57 S12-S09
     13 1324  969 (S06W19 L= 41)   213     10        0 (1305/ 948-1344/ 989)  45.0 W17-W19 S07-S04
     14  909  844 (S16E06 L= 17)   370    167        0 ( 833/ 822- 997/ 864)  11.3 E10-W00 S17-S14
     15 1838  745 (S14W64 L= 87)   447     67        1 (1815/ 713-1856/ 777)  63.4 W62-W66 S15-S11
     16 1783  645 (S20W62 L= 85)  2795    329        8 (1719/ 556-1876/ 732)  49.8 W55-W72 S26-S13
     17  233  628 (S32E66 L=316)  2836    220        2 ( 202/ 529- 268/ 731)  35.2 E83-E61 S37-S25
     18 1642  633 (S23W48 L= 70)   839    234        0 (1568/ 535-1686/ 716) -33.7 W39-W55 S28-S18
  

    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