Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1426 1692 (N30W52 L= 80)   742     82        1 (1386/1657-1473/1714) -33.7 W48-W56 N28-N32
      2  786 1625 (N48W04 L= 32)   242     71        0 ( 762/1597- 822/1653)  15.9 W00-W07 N46-N49
      3  536 1556 (N51E21 L=  7)   365     39        0 ( 514/1533- 564/1574) -11.3 E23-E18 N49-N53
      4  204 1306 (N42E62 L=326)   677    102        0 ( 156/1279- 244/1323) -11.9 E73-E54 N39-N43
      5 1209 1223 (N11W16 L= 44)   309    114        0 (1178/1175-1236/1275) -90.0 W14-W16 N07-N14
      6  509 1160 (N26E28 L=360)   394    112        2 ( 460/1133- 562/1176) -39.8 E32-E24 N22-N27
      7  171 1113 (N30E64 L=325)  1129    139        1 ( 119/1050- 221/1165) -27.9 E74-E56 N25-N34
      8  323  996 (N20E46 L=342)   472     96        1 ( 268/ 981- 371/1013) -11.3 E51-E41 N18-N21
      9 1415  506 (S37W10 L= 38)   963    165        1 (1357/ 446-1461/ 581)  37.9 W05-W14 S40-S32
     10  863  341 (S31E32 L=356)  2057    408        2 ( 731/ 222- 956/ 465)  19.5 E51-E20 S36-S24
     11 1190  371 (S39E10 L= 18)  1254    280        0 (1070/ 321-1307/ 437)   0.0 E20-W00 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