Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1224 1832 (N53E20 L=292)   552     48        0 (1193/1814-1253/1852)  -8.1 E23-E16 N50-N55
      2  649 1790 (N30E56 L=256)  1094    154        1 ( 611/1722- 705/1851) -56.3 E60-E51 N23-N35
      3  877 1757 (N33E36 L=276)  1204    242        1 ( 848/1648- 911/1834) -35.2 E44-E27 N26-N39
      4  896 1435 (N13E19 L=293)  1114    285        0 ( 865/1306- 923/1574) -60.6 E25-E15 N05-N21
      5  407 1391 (S01E49 L=263)   316     46        0 ( 382/1372- 431/1411)  45.0 E50-E47 S02-N00
      6  248 1019 (S26E54 L=258)   357     68        1 ( 227/ 987- 270/1046)  53.1 E56-E51 S28-S23
      7 1435  607 (S16W38 L=350)   375     48        1 (1414/ 584-1464/ 638) -21.8 W35-W40 S17-S14
      8 1532  468 (S20W53 L=  5)  1616    257        4 (1427/ 416-1635/ 549)  63.4 W48-W56 S26-S12
      9 1254  508 (S27W32 L=344)   934    174        2 (1171/ 489-1344/ 532)  37.9 W26-W36 S31-S23
     10  513  373 (S65E24 L=288)   559     54        1 ( 480/ 356- 544/ 394)   0.0 E29-E18 S67-S62
     11  737  303 (S61W09 L=321)  2580    296        5 ( 613/ 261- 865/ 354)  30.3 E06-W19 S68-S52
  

    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