Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  800 1866 (N46E41 L=197)   904    156        1 ( 772/1829- 821/1910) -45.0 E45-E35 N42-N51
      2 1071 1531 (N25E06 L=232)   800    168        1 (1027/1448-1127/1592)  84.8 E07-E03 N18-N29
      3  422 1436 (N09E47 L=191)   826    218        0 ( 351/1344- 503/1522)  59.7 E50-E43 N03-N15
      4 1714 1480 (N36W42 L=280)  1569    189        1 (1647/1438-1778/1517)   8.7 W35-W48 N33-N39
      5 1399 1424 (N24W16 L=254)   986    181        1 (1343/1323-1450/1508) -48.0 W11-W20 N18-N29
      6 1772 1289 (N24W46 L=284)   363     78        0 (1741/1272-1800/1309)   0.0 W43-W49 N23-N25
      7 1705 1211 (N18W41 L=278)  1103    194        1 (1674/1133-1732/1300)  79.7 W38-W42 N12-N23
      8  430 1089 (S11E39 L=199)  1518    366        1 ( 346/ 920- 504/1261)  72.3 E43-E35 S22-N00
      9 1228 1200 (N07W09 L=247)   475     79        0 (1186/1154-1265/1241) -26.6 W05-W11 N04-N09
     10  138 1198 (S07E70 L=168)   572     87        0 ( 124/1161- 152/1233) -39.8 E74-E67 S09-S04
     11 1919 1068 (N15W64 L=302)   917    149        0 (1891/1017-1941/1138)  76.0 W61-W67 N11-N19
     12  274  969 (S20E50 L=187)   302    113        2 ( 248/ 943- 307/1004) -53.1 E52-E47 S21-S17
     13  191  815 (S30E59 L=179)  1484    439        0 ( 127/ 707- 254/ 925)  47.1 E66-E51 S36-S22
     14 1373  787 (S14W25 L=263)  2752    658        4 (1061/ 725-1656/ 830)  12.7 W05-W45 S18-S08
     15  717  644 (S34E14 L=224)  1316    418        2 ( 546/ 579- 916/ 687)  24.0 E27-W00 S41-S28
     16  482  564 (S44E33 L=204)   679     88        2 ( 436/ 536- 518/ 593)  -8.1 E37-E30 S45-S41
     17 1729  475 (S23W65 L=303)  2443    312        4 (1645/ 426-1795/ 531)  26.6 W57-W71 S26-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