Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1079 1775 (N52E12 L=162)   706    136        0 (1014/1751-1144/1800)  31.0 E17-E06 N48-N55
      2  727 1568 (N29E30 L=144)  1984    392        2 ( 571/1467- 888/1644)  -9.8 E45-E16 N24-N32
      3 1606 1564 (N44W41 L=215)   812    114        1 (1545/1539-1652/1593)  15.9 W33-W47 N41-N47
      4  899 1299 (N15E12 L=162)   420    142        0 ( 846/1245- 938/1344) -35.5 E15-E08 N11-N16
      5  319 1195 (N00E50 L=124)   289     41        0 ( 300/1170- 340/1221) -14.0 E52-E47 S01-N01
      6 1882 1040 (N13W65 L=239)  4014    313        8 (1842/ 903-1932/1147) -57.3 W59-W74 N05-N19
      7 1781  941 (N06W54 L=228)   818    169        0 (1766/ 864-1796/1023) -77.5 W51-W56 N01-N10
      8  232  957 (S15E56 L=118)   666     94        2 ( 217/ 909- 251/1010) -80.5 E57-E53 S18-S12
      9  916  856 (S12E04 L=170)   379    125        0 ( 864/ 824- 962/ 897)  45.0 E06-E01 S14-S08
     10  998  678 (S21W04 L=178)  1792    514        0 ( 767/ 628-1161/ 793)  32.5 E11-W12 S27-S11
     11 1119  228 (S53W29 L=204)   876     99        0 (1085/ 193-1154/ 270) -32.5 W24-W35 S56-S49
  

    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