Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  845 1849 (N63W22 L= 39)   572     43        1 ( 816/1820- 872/1874)  26.6 W17-W27 N60-N65
      2  682 1774 (N63E02 L= 15)  1949    276        5 ( 554/1692- 792/1841)   6.5 E20-W15 N59-N66
      3  669 1481 (N40E11 L=  5)  1059    307        0 ( 542/1438- 815/1518) -26.6 E21-E01 N33-N44
      4 1636 1407 (N09W49 L= 65)  1115    194        1 (1587/1335-1682/1492)  78.7 W46-W51 N04-N14
      5 1443 1409 (N15W35 L= 52)   366     56        1 (1409/1391-1473/1431)  11.3 W32-W37 N13-N15
      6  861 1340 (N26E02 L= 15)   487     97        0 ( 808/1303- 907/1376)   0.0 E05-W01 N24-N27
      7  193 1337 (N42E63 L=314)   558     68        0 ( 160/1324- 233/1353) -10.6 E71-E55 N40-N43
      8 1832 1268 (S04W61 L= 77)   392     27        0 (1811/1246-1851/1285) -26.6 W58-W63 S04-S02
      9 1586 1178 (S01W37 L= 53)   295     49        0 (1562/1149-1611/1206)  26.6 W34-W39 S02-N00
     10  795  854 (N00E17 L=360)   214     66        0 ( 783/ 819- 807/ 887)  71.6 E17-E15 S01-N02
     11  159  821 (N11E68 L=309)  1379    147        2 ( 104/ 766- 216/ 871) -26.6 E77-E61 N07-N14
     12  790  598 (S14E24 L=352)   553    141        2 ( 763/ 532- 821/ 670)  41.2 E28-E20 S17-S10
     13 1730  564 (S44W42 L= 59)   419     85        0 (1692/ 550-1765/ 586) -21.8 W37-W47 S46-S40
     14 1540  486 (S42W21 L= 38)   844    146        1 (1462/ 467-1616/ 514) -21.8 W14-W29 S45-S38
     15  319  461 (S13E68 L=309)   754     68        0 ( 305/ 423- 335/ 501)  32.0 E73-E62 S15-S10
     16  456  460 (S15E53 L=323)   472     75        0 ( 413/ 448- 499/ 477) -15.9 E56-E49 S16-S13
     17 1315  391 (S42W01 L= 18)   532     91        0 (1276/ 364-1366/ 418)   8.1 E02-W05 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