Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  501 1403 (N41E31 L=  0)   399     81        3 ( 476/1363- 526/1428) -38.7 E33-E28 N37-N42
      2  220 1029 (N21E59 L=332)  1972    480        1 ( 159/ 853- 272/1247) -67.4 E72-E53 N10-N34
      3  626  968 (N12E25 L=  6)  1083    383        4 ( 565/ 861- 730/1125) -71.6 E28-E20 N03-N21
      4  845  975 (N08E12 L= 20)   237    104        0 ( 791/ 959- 893/ 995) -33.7 E14-E08 N06-N10
      5  382  942 (N14E44 L=347)   267     82        0 ( 362/ 905- 403/ 976) -68.2 E45-E42 N11-N16
      6 1310  573 (S27W08 L= 39)   823    258        1 (1227/ 506-1351/ 691)  45.0 W01-W12 S29-S19
      7  897  507 (S21E20 L= 11)   463    138        0 ( 859/ 457- 933/ 575)  36.9 E23-E15 S23-S17
      8  609  292 (S31E53 L=338)   446     41        0 ( 590/ 267- 625/ 318)  29.7 E56-E49 S33-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