Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  697 1732 (N63E12 L=112)   332     49        0 ( 675/1708- 721/1753)  26.6 E15-E08 N60-N64
      2  470 1332 (N35E34 L= 90)   352     87        0 ( 432/1302- 508/1359) -33.7 E38-E30 N32-N37
      3 1754  940 (S14W47 L=172)  1310    363        1 (1717/ 861-1814/1053)  65.6 W43-W53 S19-S07
      4 1890  912 (S21W64 L=189)   610     89        0 (1861/ 875-1917/ 952)  23.2 W59-W69 S22-S18
      5 1510  815 (S16W26 L=151)   241     70        0 (1495/ 781-1525/ 849)  53.1 W24-W27 S17-S13
      6 1070  650 (S16E05 L=120)   594    279        1 (1035/ 537-1160/ 712)  42.0 E09-W01 S23-S12
      7  421  538 (S13E52 L= 72)   500     85        0 ( 394/ 497- 443/ 579)  39.8 E56-E48 S15-S10
      8 1506  445 (S42W24 L=148)   313     42        0 (1488/ 422-1523/ 469)  56.3 W21-W25 S43-S39
      9 1402  369 (S45W14 L=138)   421     73        1 (1383/ 337-1431/ 403)  68.2 W11-W16 S48-S41
  

    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