Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1086 1655 (N41W27 L=257)   540     97        1 (1048/1612-1126/1690)  12.5 W22-W31 N38-N42
      2 1599 1612 (N21W60 L=290)  1580    199        5 (1491/1571-1673/1651) -13.4 W49-W70 N18-N23
      3 1541 1438 (N14W44 L=274)  1492    199        4 (1462/1358-1607/1506) -76.0 W41-W46 N07-N19
      4  825 1348 (N30E03 L=227)   788    223        1 ( 730/1295- 919/1378) -35.0 E09-W02 N24-N33
      5  740 1233 (N26E11 L=219)   248     78        0 ( 712/1193- 772/1263) -78.7 E12-E09 N22-N27
      6 1264 1100 (N04W15 L=245)  1250    350        6 (1200/ 945-1336/1246)  38.7 W07-W23 S02-N10
      7 1429 1198 (N05W27 L=258)   400    102        1 (1390/1156-1464/1244) -81.9 W26-W29 N01-N08
      8 1661  954 (S15W36 L=266)   633    163        1 (1588/ 913-1716/ 979) -16.7 W30-W40 S19-S12
      9 1866  828 (S31W55 L=285)   795     82        1 (1844/ 789-1886/ 857)  18.4 W52-W58 S33-S28
     10  799  743 (S03E20 L=210)   548    176        4 ( 766/ 676- 840/ 837)  77.5 E21-E17 S07-N02
     11  330  714 (N05E53 L=177)  1050    206        5 ( 242/ 680- 422/ 739) -17.4 E62-E45 N01-N08
     12  495  468 (S12E50 L=180)  9048      0       40 ( 358/ 314- 663/ 651)  22.4 E82-E36 S22-S02
     13  338  594 (S02E57 L=173)   362     42        1 ( 320/ 571- 360/ 616)  48.0 E59-E55 S03-N00
     14  771  395 (S23E35 L=195)  1506    362        2 ( 713/ 268- 843/ 568)  35.3 E47-E23 S30-S13
     15 1584  452 (S46W23 L=253)   501     77        1 (1565/ 421-1615/ 488) -90.0 W20-W26 S49-S42
     16  908  280 (S35E33 L=197)  1171    193        2 ( 830/ 226- 975/ 359)  19.4 E41-E24 S37-S30
     17  622  211 (S33E64 L=167)   670     79        0 ( 598/ 185- 636/ 253)  17.4 E72-E56 S35-S29
  

    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