Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  707 1615 (N38E31 L= 66)   303     93        0 ( 682/1575- 739/1647)  71.6 E32-E29 N34-N40
      2  324 1531 (N28E65 L= 32)  3630    347       10 ( 238/1432- 402/1639)  78.7 E77-E57 N20-N35
      3  401 1532 (N29E55 L= 42)   383     36        1 ( 375/1513- 421/1551)  33.7 E57-E52 N27-N30
      4 1280 1424 (N29W15 L=112)   918    279        2 (1202/1308-1333/1534) -68.2 W09-W18 N21-N36
      5 1020 1435 (N28E03 L= 94)  1281    266        2 ( 942/1324-1081/1510) -59.7 E09-W01 N20-N32
      6 1590 1298 (N23W38 L=135)   377     68        0 (1564/1264-1613/1334) -63.4 W35-W39 N20-N24
      7 1854 1138 (N14W63 L=160)   658     91        1 (1839/1098-1875/1177)  78.7 W61-W66 N11-N16
      8 1883 1081 (N10W67 L=164)   503     57        2 (1854/1063-1906/1105)   7.1 W62-W70 N09-N11
      9  177  908 (S12E65 L= 32)  2130    200        5 ( 144/ 837- 217/ 979) -41.6 E71-E60 S16-S07
     10  427  921 (S09E39 L= 58)   359    155        0 ( 364/ 872- 487/ 953)  26.6 E44-E34 S12-S06
     11 1234  663 (S19W16 L=114)  3089    749        4 ( 951/ 511-1460/ 901)  43.0 E01-W28 S31-S02
     12  772  797 (S14E14 L= 83)   306     75        1 ( 725/ 783- 812/ 817)   9.5 E17-E11 S15-S12
     13  516  745 (S19E32 L= 65)   277    135        1 ( 478/ 717- 561/ 782) -23.2 E35-E28 S20-S16
     14 1755  618 (S19W61 L=158)  3905    493        5 (1670/ 506-1830/ 762)  60.6 W53-W75 S26-S09
     15  334  621 (S30E52 L= 46)  1751    245        2 ( 253/ 531- 424/ 679)  20.0 E73-E40 S37-S24
     16 1731  509 (S27W65 L=163)   679     94        0 (1710/ 465-1758/ 544) -51.3 W62-W69 S29-S24
     17  631  424 (S41E27 L= 70)  1156    178        1 ( 555/ 378- 702/ 473)  27.9 E37-E20 S46-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