Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1143 1694 (N47E03 L=118)  1173    252        0 (1101/1618-1188/1774) -81.9 E06-W01 N41-N55
      2  619 1722 (N39E47 L= 75)   565    158        0 ( 583/1679- 645/1749) -24.0 E51-E42 N36-N40
      3  381 1656 (N30E64 L= 57)   885    175        3 ( 359/1628- 400/1695) -59.0 E66-E61 N28-N33
      4  538 1543 (N25E43 L= 78)  2094    331        3 ( 410/1472- 646/1669)  51.3 E52-E38 N19-N35
      5 1496 1580 (N44W29 L=150)   779    141        1 (1436/1553-1542/1610)  24.4 W23-W34 N41-N47
      6  575 1410 (N17E36 L= 86)   464    141        1 ( 513/1372- 617/1462)  39.8 E39-E33 N14-N20
      7 1485 1422 (N32W27 L=148)   383     98        1 (1465/1375-1501/1458)  80.5 W25-W27 N28-N34
      8  708 1251 (N09E23 L= 98)   339     84        1 ( 688/1215- 736/1288) -45.0 E24-E20 N07-N11
      9 1199 1216 (N14W08 L=129)   544    122        2 (1149/1178-1240/1271)  80.5 W05-W10 N11-N17
     10 1229 1085 (N07W11 L=133)   840    331        0 (1098/ 966-1319/1177) -29.4 W02-W18 N00-N10
     11  241  955 (S15E55 L= 66)   852    230        0 ( 205/ 863- 258/1063)  77.0 E58-E53 S21-S08
     12 1649  923 (N03W43 L=164)   639    175        3 (1626/ 869-1672/ 984) -56.3 W40-W45 N00-N06
     13  834  889 (S11E10 L=112)   343    152        0 ( 744/ 872- 901/ 910)  16.7 E15-E05 S11-S08
     14 1201  789 (S12W14 L=135)   803    188        3 (1156/ 718-1243/ 887) -68.2 W11-W17 S15-S05
     15 1627  734 (S09W45 L=166)   308    118        0 (1599/ 665-1658/ 782)  76.0 W43-W46 S13-S05
     16 1735  693 (S10W56 L=177)   494     80        1 (1713/ 650-1755/ 733)  90.0 W54-W57 S12-S07
     17  694  561 (S34E16 L=105)   293     51        0 ( 668/ 543- 724/ 578)   0.0 E18-E13 S35-S32
  

    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