Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  461 1327 (N34E31 L=178)  1012    180        3 ( 387/1281- 525/1367) -28.6 E36-E25 N29-N35
      2  549  736 (S03E35 L=174)  2661    694        8 ( 455/ 528- 700/ 836)  29.5 E48-E25 S13-N02
      3 1063  771 (S12E04 L=205)   208     81        1 (1013/ 759-1103/ 785) -21.8 E06-E01 S13-S10
      4  293  723 (N01E56 L=153)   341     69        1 ( 276/ 695- 313/ 759)  71.6 E57-E53 N00-N03
      5  700  579 (S15E30 L=178)  2175    488        5 ( 590/ 439- 771/ 725)  35.0 E42-E22 S21-S07
      6 1051  675 (S17E07 L=202)   253     57        1 (1019/ 649-1071/ 703)  33.7 E08-E04 S18-S15
      7  833  461 (S25E26 L=182)   918    284        6 ( 755/ 368- 921/ 518)  19.4 E35-E18 S29-S22
      8  472  467 (S17E51 L=158)   333     60        0 ( 444/ 436- 492/ 494) -90.0 E52-E49 S19-S14
      9  495  392 (S22E53 L=155)   367     68        0 ( 464/ 376- 528/ 410) -11.3 E55-E50 S23-S20
  

    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