Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1247 1784 (N56E03 L= 35)  1738    228        2 (1122/1747-1338/1813)  28.8 E14-W06 N50-N63
      2  564 1481 (N17E41 L=358)   268     83        2 ( 522/1457- 592/1499)  21.8 E43-E37 N15-N18
      3 1181 1391 (N24W02 L= 41)   255     70        0 (1143/1372-1212/1412) -14.0 W00-W04 N22-N24
      4 1827 1079 (N18W57 L= 95)  3833    483        7 (1750/ 935-1914/1213) -58.0 W48-W71 N10-N26
      5 1724 1071 (N16W46 L= 85)   816    168        2 (1691/ 994-1750/1139) -58.0 W42-W48 N10-N19
      6  329 1062 (S12E47 L=351)   371    100        0 ( 297/1016- 368/1098)  50.2 E49-E44 S15-S09
      7  931  766 (S19E00 L= 38)   424    159        0 ( 865/ 742- 978/ 809)  26.6 E04-W02 S20-S15
      8  188  747 (S35E68 L=331)   928     79        2 ( 164/ 711- 216/ 783) -38.7 E72-E62 S37-S32
      9 1380  557 (S22W34 L= 72)   746    191        1 (1312/ 515-1451/ 636)  68.2 W30-W36 S25-S15
     10 1510  374 (S31W56 L= 95)  2301    303        4 (1363/ 318-1631/ 420)  28.8 W46-W67 S37-S24
     11 1310  201 (S48W59 L= 98)   583     47        1 (1283/ 179-1335/ 229)   5.7 W53-W66 S49-S46
  

    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