Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  628 1639 (N30E36 L=192)   402     67        1 ( 595/1617- 662/1656)  -9.5 E38-E32 N28-N30
      2  778 1520 (N23E22 L=206)   273     76        1 ( 749/1492- 816/1550) -21.8 E23-E18 N20-N24
      3  461 1244 (N03E38 L=190)   381     81        0 ( 444/1196- 483/1291)  90.0 E38-E36 N00-N05
      4  823  921 (S14E11 L=217)   224     69        1 ( 785/ 905- 859/ 942)  26.6 E13-E09 S14-S11
      5 1382  792 (S15W25 L=253)  1227    253        4 (1272/ 719-1497/ 869)  45.0 W18-W31 S20-S09
      6  647  777 (S24E22 L=206)   203     63        0 ( 613/ 764- 685/ 787)  11.3 E24-E19 S24-S22
      7  802  689 (S28E11 L=217)   368    130        0 ( 738/ 647- 855/ 732)  45.0 E15-E07 S31-S24
      8 1576  651 (S21W42 L=270)   403    110        0 (1529/ 600-1604/ 706)  66.8 W39-W43 S24-S17
      9  652  343 (S55E26 L=201)   514    102        1 ( 626/ 311- 690/ 371)  24.0 E30-E21 S57-S52
  

    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