Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  338 1112 (S15E46 L=207)  1110    246        3 ( 272/1024- 404/1214)  65.2 E50-E43 S21-S08
      2  240 1143 (S16E56 L=197)   365     59        0 ( 223/1103- 260/1172)  90.0 E57-E54 S17-S13
      3  174  820 (S37E59 L=195)  1354    305        0 ( 128/ 786- 210/ 865)  31.0 E64-E54 S39-S33
      4  389  789 (S34E37 L=217)   708    161        2 ( 320/ 750- 447/ 833)  -9.5 E43-E31 S35-S32
      5  787  772 (S24E07 L=246)  1689    469        4 ( 547/ 695- 974/ 822)   9.8 E24-W05 S28-S21
      6 1349  677 (S15W29 L=283)  1697    484        6 (1146/ 637-1575/ 733)  17.8 W15-W43 S18-S07
      7  536  554 (S46E22 L=232)   568    100        2 ( 509/ 511- 573/ 596) -26.6 E24-E18 S48-S42
      8  667  528 (S44E10 L=244)   299     97        0 ( 622/ 506- 704/ 551)  45.0 E13-E07 S46-S40
      9 1404  289 (S35W54 L=307)   523     94        1 (1361/ 273-1453/ 309)  51.3 W51-W55 S37-S32
  

    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