Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  929 1741 (N51W20 L=143)   532     89        1 ( 904/1695- 961/1779)  29.1 W15-W24 N47-N53
      2 1292 1713 (N36W43 L=165)   522    120        1 (1238/1678-1351/1738) -11.3 W38-W48 N35-N37
      3 1251 1599 (N31W33 L=155)   248     34        0 (1221/1587-1278/1614) -14.0 W30-W34 N29-N31
      4  809 1356 (N29E04 L=118)   272     48        0 ( 774/1338- 838/1376) -36.9 E06-E02 N27-N30
      5 1734 1275 (S01W51 L=174)   941    137        1 (1691/1238-1798/1317) -45.0 W46-W56 S04-N02
      6  345 1077 (N24E43 L= 80)  2241    415        4 ( 192/ 989- 518/1209) -33.2 E59-E30 N15-N34
      7 1811  947 (S21W50 L=172)   528    122        2 (1783/ 891-1830/1009)  50.2 W46-W52 S23-S17
      8 1918  941 (S26W64 L=186)   652     54        0 (1888/ 923-1949/ 959) -10.3 W59-W70 S26-S23
      9  248  808 (N10E57 L= 65)  1131    149        4 ( 206/ 741- 295/ 844)  71.6 E62-E52 N06-N12
     10 1269  732 (S17W06 L=129)   724    320        0 (1165/ 643-1393/ 833)  14.0 E01-W15 S20-S14
     11  905  770 (S06E13 L=109)  1211    313        2 ( 840/ 665-1032/ 816)   4.4 E19-E06 S10-S02
     12 1595  651 (S32W27 L=149)  1400    333        0 (1499/ 556-1709/ 785)  21.8 W17-W37 S36-S26
     13  279  710 (N04E57 L= 65)  1840    236        6 ( 198/ 629- 352/ 781)   0.0 E68-E50 N00-N07
     14  475  591 (S06E45 L= 77)   343     57        0 ( 449/ 574- 506/ 611)   0.0 E47-E42 S07-S05
     15  588  503 (S14E41 L= 81)   871    205        1 ( 523/ 436- 640/ 591)  26.6 E49-E34 S16-S09
     16  945  237 (S41E33 L= 89)   400     66        1 ( 925/ 205- 961/ 280)  36.9 E37-E29 S43-S37
     17  838  162 (S44E50 L= 72)   564     74        0 ( 803/ 142- 871/ 190)   5.2 E55-E43 S45-S42
  

    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