Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1766 (N49W47 L=328)   579     83        0 (1269/1735-1349/1783) -14.0 W42-W53 N46-N49
      2 1428 1695 (N40W52 L=332)  1090    126        3 (1384/1655-1475/1740)  81.9 W46-W55 N36-N43
      3 1386 1600 (N35W40 L=321)   567    170        1 (1330/1532-1421/1636)  82.9 W36-W43 N30-N38
      4 1326 1563 (N34W34 L=314)   281    116        1 (1294/1524-1363/1611) -81.9 W32-W35 N30-N37
      5  905 1185 (N18E05 L=275)   656    316        0 ( 771/1139-1046/1229) -31.0 E13-W03 N13-N22
      6  403 1150 (N22E41 L=239)   591    181        2 ( 354/1065- 432/1212) -63.4 E46-E38 N16-N26
      7 1829  952 (S14W57 L=337)  2979    527        5 (1760/ 729-1882/1147)  63.4 W48-W66 S27-S02
      8  158 1118 (N22E70 L=210)   521     51        0 ( 139/1100- 186/1141) -20.6 E73-E65 N20-N23
      9  784 1015 (N09E15 L=265)   388    119        2 ( 734/ 979- 833/1050) -50.2 E17-E12 N06-N12
     10 1431  930 (S06W23 L=304)   482    169        1 (1394/ 847-1465/ 995)  53.1 W19-W26 S10-S02
     11  543  469 (S22E44 L=236)   513    179        0 ( 522/ 411- 562/ 565)  57.5 E49-E40 S26-S15
     12  408  454 (S22E59 L=222)   429     79        1 ( 386/ 423- 430/ 494)  29.1 E64-E54 S24-S19
     13 1070  196 (S55E18 L=262)   439     63        1 (1042/ 180-1104/ 219)   0.0 E21-E13 S56-S51
  

    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