Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1328 1812 (N55W06 L=193)   707    108        1 (1285/1794-1366/1837)  32.0 W01-W09 N51-N59
      2 1331 1731 (N46W07 L=195)   495     83        1 (1304/1702-1351/1766)  90.0 W04-W08 N44-N50
      3  473 1469 (N11E44 L=144)   415    104        0 ( 453/1428- 494/1524) -71.6 E45-E41 N08-N14
      4 1134 1341 (N14W01 L=188)   331    104        0 (1103/1291-1154/1396)  90.0 W00-W01 N10-N17
      5 1768 1205 (N20W46 L=234)  1458    318        1 (1734/1035-1814/1351) -67.2 W42-W52 N10-N29
      6 1239  953 (S07W14 L=202)   963    214        2 (1183/ 862-1269/1068) -80.5 W11-W15 S12-N00
      7  599  851 (S24E24 L=164)   668    220        1 ( 491/ 793- 686/ 913)  39.8 E31-E18 S29-S19
      8  288  884 (S26E49 L=139)   337     65        0 ( 263/ 860- 312/ 907)  36.9 E50-E46 S27-S24
      9 1093  769 (S20W09 L=197)   531    155        1 (1020/ 721-1170/ 828)  56.3 W05-W12 S24-S15
     10  878  760 (S25E04 L=184)   432    105        1 ( 830/ 718- 910/ 806)  66.8 E06-E01 S28-S21
     11  628  619 (S38E20 L=168)   832    230        0 ( 517/ 602- 723/ 643)  16.4 E29-E12 S40-S34
     12 1519  499 (S26W46 L=234)  4955    942       10 (1130/ 424-1844/ 588)  24.5 W21-W78 S39-S13
  

    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