Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  600 1675 (N25E48 L=247)  1663    455        4 ( 534/1552- 638/1784) -35.0 E60-E39 N18-N32
      2  818 1734 (N35E36 L=259)   595    136        0 ( 771/1704- 868/1760)  -5.7 E41-E31 N33-N36
      3  269 1486 (N07E70 L=226)   622    102        1 ( 255/1449- 287/1516) -38.7 E73-E66 N04-N08
      4  551 1404 (N08E39 L=256)   356     56        2 ( 529/1377- 578/1439) -76.0 E40-E36 N05-N09
      5  200 1239 (S09E65 L=230)   562    100        1 ( 173/1206- 219/1267) -20.6 E70-E62 S10-S07
      6  215 1156 (S13E60 L=235)  2120    337        5 ( 157/1062- 272/1267)  70.3 E65-E54 S19-S05
      7  652  998 (S13E21 L=274)   432     73        2 ( 617/ 963- 689/1037)  76.0 E23-E18 S15-S11
      8  191  837 (S34E57 L=238)   528     81        1 ( 168/ 811- 215/ 863) -26.6 E60-E54 S35-S32
      9  621  755 (S29E19 L=276)   312     61        0 ( 592/ 723- 649/ 779)  59.0 E20-E17 S31-S26
     10  471  620 (S42E30 L=266)   687    115        1 ( 429/ 584- 512/ 662)  45.0 E33-E25 S45-S38
     11 1556  367 (S27W59 L=354)   450    103        1 (1525/ 352-1589/ 388)   8.1 W55-W62 S27-S25
  

    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