Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1268 1582 (N40W15 L=299)  7352   1244       12 ( 731/1435-1769/1727)  -6.4 E33-W65 N30-N47
      2  911 1192 (N12E08 L=275)  1871    422        2 ( 880/ 992- 934/1382) -77.7 E11-E05 N00-N23
      3  275 1181 (N07E55 L=229)  2423    392        3 ( 229/1012- 311/1362) -73.3 E62-E49 S03-N17
      4 1663 1276 (N21W44 L=327)  1064    236        1 (1605/1221-1716/1331) -45.0 W39-W49 N17-N24
      5  503 1284 (N15E37 L=247)   341    100        1 ( 461/1256- 556/1309)  15.9 E39-E32 N12-N16
      6  418 1207 (N09E42 L=241)   651    146        2 ( 392/1143- 458/1274) -69.4 E44-E37 N05-N13
      7  476  765 (S17E36 L=248)  1494    425        3 ( 431/ 570- 509/ 949)  72.3 E40-E32 S30-S05
      8  190  892 (S12E63 L=220)   972     93        1 ( 151/ 846- 223/ 926) -26.6 E68-E58 S14-S09
      9  909  624 (S24E05 L=278)   469    130        0 ( 847/ 596- 970/ 660)  26.6 E09-E01 S25-S20
     10 1760  584 (S22W63 L=346)  2594    217        4 (1663/ 534-1856/ 622)   6.3 W51-W79 S25-S19
     11 1471  447 (S33W39 L=322)  2345    393        4 (1270/ 371-1617/ 503)  22.9 W24-W51 S39-S28
     12 1083  318 (S46W11 L=294)   394     64        0 (1046/ 301-1116/ 334)  18.4 W07-W13 S47-S43
  

    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