Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1034 1779 (N56E06 L=316)   451     82        1 ( 989/1768-1083/1792)  11.3 E11-E01 N54-N58
      2 1256 1581 (N40W15 L=337)   418     93        1 (1202/1561-1290/1602)  15.9 W10-W17 N38-N42
      3  520 1407 (N23E39 L=283)   321     72        1 ( 482/1390- 559/1428)  -8.1 E42-E35 N22-N24
      4 1169 1277 (N19W08 L=330)   647    268        1 (1064/1195-1275/1330) -25.0 W00-W15 N14-N21
      5  720 1232 (N13E21 L=301)   578    142        0 ( 657/1201- 776/1278) -14.0 E25-E17 N11-N15
      6 1658 1040 (N06W43 L=  5)   485     47        2 (1637/1007-1678/1071)  53.1 W41-W44 N04-N08
      7  691  723 (S18E20 L=302)   465    140        0 ( 656/ 660- 724/ 786)  63.4 E22-E17 S22-S14
      8 1620  706 (S15W44 L=  6)   386     62        1 (1595/ 673-1643/ 741)  76.0 W42-W45 S16-S12
      9 1711  651 (S18W54 L= 16)   347     34        0 (1694/ 632-1729/ 676) -45.0 W52-W56 S19-S16
     10  879  641 (S23E07 L=315)   340    136        0 ( 822/ 627- 937/ 669)  -7.1 E11-E03 S23-S21
     11  345  522 (S36E57 L=265)   949    127        1 ( 290/ 503- 408/ 549)   0.0 E65-E48 S37-S33
  

    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