Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1403 1133 (N03W23 L= 57)   257     50        0 (1373/1105-1436/1157) -36.9 W21-W25 N01-N04
      2  273  869 (S10E53 L=341)  1075    190        1 ( 254/ 771- 302/ 968)  80.5 E54-E50 S16-S04
      3 1895  902 (S10W66 L=100)   549    123        0 (1857/ 884-1927/ 920)   5.7 W60-W70 S10-S08
      4  591  848 (S13E28 L=  6)   233     82        0 ( 551/ 826- 637/ 863)   0.0 E30-E24 S14-S11
      5  135  814 (S12E70 L=324)  1016     98        0 ( 109/ 783- 160/ 847) -29.7 E74-E66 S14-S10
      6  581  661 (S24E31 L=  3)   400     43        1 ( 555/ 631- 613/ 688) -33.7 E32-E28 S26-S22
      7  362  566 (S30E51 L=343)  1483    261        2 ( 310/ 498- 408/ 640)  32.7 E58-E44 S33-S24
      8 1673  521 (S34W50 L= 84)   432     64        0 (1637/ 500-1698/ 549)  45.0 W47-W52 S35-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