Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  721 1792 (N63W03 L= 21)  2493    344        7 ( 587/1706- 864/1873)   5.6 E15-W27 N59-N65
      2  657 1473 (N40E12 L=  5)  1151    325        1 ( 508/1429- 801/1512) -25.6 E25-E02 N33-N44
      3 1627 1399 (N09W47 L= 65)   929    187        2 (1572/1328-1672/1487)  84.8 W45-W50 N04-N15
      4 1428 1400 (N15W34 L= 52)   363     65        1 (1389/1384-1468/1422) -21.8 W31-W36 N13-N15
      5  842 1332 (N26E03 L= 15)   489    121        0 ( 800/1285- 894/1369)  18.4 E06-W00 N23-N27
      6  190 1334 (N42E64 L=314)   599     78        0 ( 155/1318- 229/1352) -10.0 E73-E56 N40-N43
      7 1831 1267 (S04W60 L= 78)   435     30        1 (1813/1242-1851/1290)  33.7 W58-W63 S04-S02
      8 1572 1173 (S01W36 L= 53)   311     42        0 (1548/1145-1596/1202)  14.0 W33-W37 S02-N00
      9  857 1043 (N09E09 L=  9)   204     80        1 ( 816/1023- 896/1063) -53.1 E10-E06 N07-N11
     10  152  819 (N11E69 L=309)  1401    157        4 ( 102/ 770- 210/ 868) -25.0 E77-E62 N07-N14
     11  756  554 (S16E28 L=350)   812    232        3 ( 704/ 448- 803/ 656)  40.2 E35-E21 S21-S10
     12 1721  560 (S44W41 L= 59)   470     87        1 (1682/ 543-1762/ 583) -10.3 W35-W46 S45-S40
     13 1524  482 (S42W20 L= 38)   822    156        2 (1441/ 464-1600/ 509) -18.4 W12-W27 S45-S39
     14  311  465 (S12E69 L=309)   779     71        1 ( 287/ 432- 329/ 502)  26.6 E73-E63 S14-S09
     15  451  456 (S15E54 L=324)   421     78        1 ( 411/ 438- 494/ 468) -18.4 E57-E50 S16-S13
     16 1307  387 (S42W00 L= 18)   488     82        0 (1267/ 363-1358/ 413)   0.0 E02-W04 S43-S40
  

    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