Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  608 1819 (N42E56 L= 92)  1313      0       -1 ( 584/1781- 632/1859) -22.8 E68-E49 N38-N46
      2  959 1676 (N34E16 L=132)   324     61        2 ( 926/1652- 990/1694)  53.1 E17-E13 N32-N36
      3  373 1638 (N25E64 L= 84)   709     66        2 ( 356/1599- 391/1676) -26.6 E71-E59 N22-N28
      4  911 1385 (N14E12 L=136)   396     95        2 ( 876/1356- 948/1419)  59.0 E14-E09 N11-N16
      5 1784 1308 (N25W51 L=199)  1425    189        0 (1731/1212-1833/1378) -45.0 W45-W55 N19-N28
      6 1514 1313 (N19W27 L=175)   285     44        0 (1493/1285-1531/1345)  76.0 W25-W27 N17-N21
      7 1081 1253 (N08W00 L=148)   199     58        0 (1055/1220-1102/1276) -56.3 E01-W01 N06-N09
      8  304 1237 (S02E52 L= 96)   305     64        0 ( 293/1196- 318/1265)  90.0 E52-E50 S04-N00
      9 1842 1174 (N17W56 L=204)   908     92        1 (1797/1137-1878/1208)   0.0 W51-W60 N15-N20
     10 1274 1101 (N02W14 L=162)   213     89        2 (1236/1068-1303/1120) -14.0 W11-W15 N00-N03
     11 1364 1065 (N01W20 L=168)   238     96        0 (1322/1041-1400/1106)  45.0 W17-W21 N00-N04
     12 1702 1011 (N04W44 L=192)   337     99        0 (1669/ 958-1730/1049) -33.7 W40-W46 N01-N06
     13 1312  963 (S05W18 L=166)   501    224        0 (1208/ 915-1408/ 991)  -4.4 W11-W24 S06-S04
     14 1340  786 (S15W23 L=171)   688    255        1 (1281/ 681-1406/ 880)  77.9 W21-W25 S22-S08
     15 1584  700 (S16W42 L=190)   260     53        0 (1566/ 685-1611/ 728) -14.0 W40-W44 S16-S14
     16 1385  637 (S23W30 L=178)   347     80        1 (1342/ 621-1425/ 661)   0.0 W26-W33 S23-S21
     17 1807  578 (S17W71 L=220)   924     83        0 (1762/ 558-1851/ 596)  15.9 W66-W80 S19-S15
     18 1607  446 (S30W58 L=206)   394     46        1 (1580/ 427-1632/ 465)  21.8 W54-W60 S30-S28
  

    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