Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1171 1802 (N52W44 L= 61)   349     37        1 (1150/1779-1197/1820)   0.0 W40-W48 N50-N53
      2  244 1267 (N36E59 L=318)   638     93        0 ( 214/1220- 277/1304) -28.6 E64-E53 N32-N38
      3  803 1074 (N15E12 L=  4)   944    268        3 ( 769/ 993- 860/1208) -90.0 E13-E09 N08-N22
      4  393 1073 (N22E41 L=335)   639    256        1 ( 305/1044- 494/1136) -26.6 E49-E33 N18-N26
      5  940  948 (N05E07 L= 10)   200     50        0 ( 908/ 934- 967/ 963) -33.7 E08-E04 N03-N05
      6  442  880 (N10E40 L=336)   274     64        2 ( 404/ 856- 469/ 903) -36.9 E42-E38 N07-N11
      7 1541  866 (S14W28 L= 44)   232     74        0 (1520/ 834-1559/ 902)  45.0 W26-W29 S15-S12
      8  217  816 (N08E63 L=313)   807    137        2 ( 179/ 792- 279/ 866) -21.8 E67-E56 N06-N11
      9 1114  584 (S21E04 L= 12)   532    159        0 (1070/ 520-1154/ 653)  41.2 E08-W00 S24-S17
     10  743  328 (S31E39 L=338)   434     90        0 ( 726/ 284- 767/ 355)  56.3 E41-E35 S34-S28
     11  595  274 (S33E57 L=319)  1183    125        3 ( 552/ 217- 625/ 323)  57.5 E63-E52 S38-S27
  

    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