Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  988 1907 (N63W50 L=199)   521     61        0 ( 971/1875-1006/1941)  12.8 W41-W63 N60-N65
      2  429 1462 (N47E34 L=116)   423     49        2 ( 408/1433- 456/1488)  45.0 E36-E30 N44-N48
      3 1015 1326 (N22W08 L=157)   215     80        0 ( 975/1303-1055/1339) -11.3 W04-W09 N20-N22
      4  341 1116 (N26E43 L=106)   660    204        0 ( 301/1043- 420/1161)  15.3 E47-E36 N22-N28
      5 1700  914 (S20W40 L=189)   371    140        0 (1661/ 868-1756/ 957)  11.3 W35-W45 S21-S18
      6 1623  683 (S32W30 L=180)   482    130        0 (1566/ 647-1679/ 727)   5.2 W24-W35 S32-S30
      7  839  564 (S17E23 L=126)   257     94        0 ( 794/ 538- 875/ 595)   9.5 E26-E20 S17-S15
      8  437  526 (S10E52 L= 97)   815    151        0 ( 414/ 445- 464/ 583)  56.3 E56-E47 S15-S06
      9 1410  491 (S38W10 L=160)   654    144        0 (1347/ 451-1475/ 537)  10.3 W04-W15 S39-S36
     10 1236  406 (S38E05 L=145)   247     33        2 (1212/ 393-1263/ 424)   0.0 E06-E02 S39-S36
     11  886  303 (S35E33 L=117)   290     60        0 ( 866/ 276- 902/ 336)  29.7 E36-E29 S36-S32
     12  558  284 (S28E61 L= 89)  1479    148        3 ( 505/ 247- 618/ 334) -90.0 E65-E54 S31-S23
  

    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