Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  574 1665 (N46E45 L=342)   483     91        1 ( 542/1630- 600/1692) -20.0 E50-E39 N42-N47
      2  627 1588 (N40E34 L=353)   295     61        1 ( 606/1563- 654/1609)  63.4 E35-E32 N37-N41
      3 1121 1387 (N26W06 L= 33)   261     96        1 (1095/1333-1140/1418) -71.6 W04-W07 N22-N28
      4  289 1278 (N17E56 L=331)  1042    214        1 ( 258/1191- 330/1372) -42.5 E62-E49 N12-N23
      5  748 1206 (N14E18 L=  9)   783    276        0 ( 692/1100- 827/1306) -55.3 E22-E12 N07-N20
      6  784  878 (S06E15 L= 12)   238     70        0 ( 772/ 837- 797/ 915)  78.7 E15-E14 S08-S03
      7  540  783 (S12E32 L=355)   274    105        1 ( 515/ 728- 554/ 823) -71.6 E33-E30 S15-S09
      8 1689  500 (S32W58 L= 85)   912     82        2 (1656/ 465-1711/ 539) -68.2 W53-W62 S34-S29
  

    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