Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  232 1087 (S14E56 L=293)  2258    382       10 ( 151/1004- 285/1176)  79.7 E64-E50 S19-S08
      2  436 1079 (S12E38 L=311)   238     78        0 ( 410/1043- 470/1113)  53.1 E39-E35 S14-S10
      3  633 1058 (S11E24 L=325)   226     97        0 ( 596/1009- 667/1097) -38.7 E26-E21 S13-S09
      4 1814 1042 (N11W52 L= 41)   363     70        0 (1795/1010-1832/1077) -36.9 W50-W54 N09-N12
      5 1192  898 (S11W12 L=  1)   232     87        0 (1148/ 863-1225/ 943) -26.6 W08-W14 S12-S08
      6 1891  877 (N04W66 L= 55)  1641    276        3 (1840/ 794-1927/ 927) -14.0 W58-W75 N00-N06
      7 1125  775 (S19W11 L=360)   482    128        1 (1062/ 755-1202/ 807)  35.5 W06-W14 S21-S15
      8 1501  774 (S12W34 L= 23)   569    176        2 (1433/ 727-1565/ 803)   5.7 W29-W39 S13-S10
      9 1610  736 (S11W43 L= 32)   279     65        0 (1594/ 705-1624/ 773) -45.0 W40-W44 S12-S09
     10  335  632 (S42E46 L=303)   930    291        1 ( 288/ 550- 390/ 717)  42.7 E52-E39 S47-S35
     11 1114  653 (S27W13 L=  2)   648    180        0 (1037/ 611-1202/ 710)  48.4 W08-W17 S30-S21
     12 1621  641 (S16W47 L= 36)   499    147        3 (1566/ 609-1680/ 682)  32.0 W43-W51 S18-S12
     13 1828  621 (S11W69 L= 58)   522    130        0 (1792/ 604-1859/ 640)  24.0 W64-W73 S13-S09
     14 1810  556 (S15W71 L= 60)   990    177        4 (1778/ 497-1837/ 596)  74.1 W67-W73 S19-S12
     15 1596  531 (S23W50 L= 39)  1078    291        2 (1492/ 469-1671/ 578)  10.0 W41-W58 S25-S20
     16 1259  509 (S32W27 L= 16)  1414    254        6 (1127/ 475-1403/ 546)  24.0 W18-W36 S35-S27
  

    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