Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1771 1285 (N16W58 L=336)   489     93        1 (1728/1254-1802/1308) -29.7 W53-W60 N13-N17
      2 1525 1202 (N13W34 L=312)   487    132        0 (1467/1163-1570/1248) -45.0 W30-W37 N10-N16
      3  682  856 (S05E22 L=256)   413    140        0 ( 606/ 822- 732/ 879) -20.6 E27-E19 S07-S03
      4  277  840 (S06E55 L=224)   622    127        1 ( 226/ 821- 330/ 878) -16.7 E59-E49 S07-S04
      5  727  606 (S22E22 L=257)   339    173        0 ( 690/ 551- 757/ 661) -82.9 E23-E19 S25-S17
      6  321  604 (S22E57 L=222)   433    142        2 ( 288/ 572- 375/ 642) -26.6 E59-E51 S24-S19
      7  610  600 (S22E30 L=248)   461    156        0 ( 561/ 573- 660/ 638)   8.1 E33-E26 S23-S19
      8  917  475 (S31E10 L=269)   330     79        1 ( 879/ 448- 960/ 507)  26.6 E12-E06 S33-S28
      9 1327  421 (S38W21 L=300)   765    154        1 (1294/ 353-1372/ 482)  73.3 W19-W24 S43-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