Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1219 1634 (N37E07 L= 16)  2021    451        4 (1141/1473-1305/1799) -39.8 E21-W03 N28-N48
      2  414 1684 (N19E66 L=317)  2556    290        4 ( 363/1635- 468/1732)  71.6 E70-E62 N14-N23
      3  608 1580 (N16E42 L=340)   587    167        0 ( 585/1523- 641/1649) -45.0 E46-E38 N13-N20
      4  502 1591 (N14E51 L=332)   386     63        0 ( 468/1574- 531/1606)  11.3 E53-E48 N12-N15
      5  369 1513 (N07E59 L=324)   731    114        1 ( 321/1456- 407/1548)  66.8 E61-E56 N02-N09
      6  458 1493 (N07E50 L=333)   481     70        1 ( 412/1450- 487/1526)  59.0 E51-E47 N04-N09
      7 1681 1235 (N26W36 L= 58)   665    151        0 (1654/1161-1699/1304) -58.0 W32-W37 N21-N29
      8  840 1047 (S09E11 L= 12)   200     29        0 ( 825/1024- 857/1074) -33.7 E12-E09 S09-S07
      9 1837  909 (N13W57 L= 80)  2478    389        6 (1798/ 763-1882/1042) -39.1 W50-W68 N06-N20
     10 1524  937 (N04W31 L= 54)   235     15        0 (1505/ 917-1542/ 958)  71.6 W30-W32 N02-N05
     11  189  798 (S39E57 L=326)   715    184        1 ( 144/ 784- 220/ 820)  18.4 E62-E53 S40-S37
     12 1269  782 (S12W21 L= 43)   536    189        0 (1203/ 712-1363/ 817)  63.4 W18-W25 S17-S07
     13  995  800 (S18W05 L= 27)   269     93        0 ( 940/ 787-1042/ 811)  26.6 W01-W07 S19-S16
     14 1017  655 (S26W10 L= 33)   677    165        1 ( 925/ 637-1101/ 678)  31.0 W04-W15 S28-S22
     15 1185  507 (S29W26 L= 49)   332     66        0 (1147/ 487-1225/ 527)   9.5 W23-W29 S30-S28
     16  916  465 (S41W12 L= 34)   283     52        0 ( 897/ 444- 943/ 493) -26.6 W10-W14 S42-S38
     17 1007  460 (S38W18 L= 41)   249     36        0 ( 980/ 445-1035/ 475)  26.6 W15-W19 S39-S36
  

    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