Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  822 1615 (N49W00 L= 98)   825    137        3 ( 771/1559- 862/1650) -23.2 E04-W03 N44-N51
      2  839 1329 (N29E05 L= 92)   413    147        0 ( 806/1275- 893/1382)  32.0 E08-W00 N26-N31
      3  435 1146 (N25E38 L= 60)  1564    539        1 ( 288/1057- 645/1220) -22.5 E51-E22 N16-N30
      4  223 1065 (N22E59 L= 38)   792    125        2 ( 191/1000- 254/1124)  45.0 E63-E55 N18-N25
      5  699  657 (S09E28 L= 70)   257     37        1 ( 679/ 636- 719/ 681) -71.6 E29-E26 S10-S07
      6 1024  452 (S28E13 L= 85)   576    144        1 ( 974/ 399-1074/ 506)  24.0 E17-E08 S31-S25
      7  514  427 (S21E51 L= 47)   629    107        1 ( 485/ 395- 559/ 456) -59.0 E54-E48 S23-S18
      8  641  365 (S27E44 L= 54)  1248    234        1 ( 591/ 277- 718/ 447)  82.4 E46-E40 S35-S20
      9  738  244 (S39E46 L= 52)   410     54        0 ( 720/ 219- 753/ 275)  26.6 E50-E42 S41-S36
  

    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