Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1156 1569 (N31W09 L=  2)   672    131        1 (1084/1516-1221/1616) -33.7 W04-W13 N27-N34
      2  404 1365 (N18E42 L=311)   547    114        0 ( 384/1316- 423/1418) -60.3 E44-E40 N14-N21
      3 1830 1324 (N17W59 L= 52)   531     95        0 (1791/1307-1859/1338)   0.0 W55-W62 N15-N17
      4  268 1232 (N10E53 L=300)   533    111        1 ( 228/1218- 311/1249) -14.0 E56-E48 N09-N11
      5  545 1188 (N07E30 L=323)   273     76        0 ( 520/1140- 568/1231)  63.4 E31-E28 N03-N09
      6  857 1090 (N00E10 L=343)   306     65        1 ( 815/1061- 880/1118)   0.0 E12-E08 S01-N02
      7  846  911 (S10E11 L=342)   543    165        0 ( 818/ 820- 883/ 998)  68.2 E12-E08 S15-S05
      8 1749  801 (S15W51 L= 44)   462    147        0 (1693/ 785-1801/ 818)   0.0 W45-W55 S16-S14
      9  296  667 (S24E54 L=299)  1475    269        1 ( 250/ 593- 322/ 772)  54.0 E60-E49 S28-S17
     10  564  658 (S25E31 L=322)   548     91        1 ( 543/ 613- 588/ 710)  63.4 E33-E29 S28-S22
     11 1149  633 (S27W08 L=  2)   418     95        0 (1122/ 583-1167/ 682) -71.6 W06-W09 S30-S24
     12  505  487 (S37E40 L=313)   359     57        1 ( 480/ 464- 531/ 513) -26.6 E42-E38 S38-S34
  

    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