Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  536 1512 (N34E40 L=162)   844    152        1 ( 463/1486- 613/1532)  -3.8 E47-E32 N31-N35
      2  701 1474 (N32E25 L=177)   325     74        0 ( 667/1451- 736/1508) -23.2 E28-E21 N30-N33
      3  746  687 (S18E18 L=184)   439    206        0 ( 654/ 664- 838/ 719)  -8.7 E24-E11 S19-S16
      4  402  579 (S27E48 L=154)  2332    332        3 ( 319/ 422- 476/ 688)  27.7 E78-E38 S40-S19
      5 1708  650 (S21W53 L=254)   422     37        0 (1688/ 629-1729/ 674) -31.0 W50-W55 S22-S19
      6  808  612 (S23E14 L=187)   392    184        0 ( 759/ 546- 883/ 653)  35.0 E18-E08 S27-S20
      7 1576  414 (S38W50 L=252)   684    102        0 (1549/ 366-1598/ 455) -63.4 W46-W52 S42-S34
      8 1463  363 (S42W41 L=242)  1104    112        3 (1408/ 339-1522/ 396)   4.8 W35-W47 S44-S39
      9 1231  308 (S47W20 L=222)  1063    197        3 (1145/ 287-1320/ 344)   7.1 W12-W28 S48-S43
  

    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