Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1395 1743 (N38W50 L= 67)   513     53        1 (1374/1702-1412/1772)  26.6 W45-W54 N35-N39
      2 1555 1714 (N32W63 L= 80)  1301    131        0 (1497/1688-1599/1749)  -8.7 W55-W70 N29-N34
      3  600 1697 (N56E21 L=356)  1462    116        2 ( 558/1637- 653/1742)  71.6 E27-E14 N50-N60
      4  909 1596 (N39W05 L= 23)  1604    273        2 ( 795/1530-1043/1647) -42.9 E02-W12 N31-N45
      5 1308 1596 (N30W33 L= 50)   308     62        0 (1271/1574-1338/1614) -45.0 W30-W34 N28-N31
      6 1034 1390 (N23W08 L= 25)   451     70        2 ( 998/1361-1072/1412)   0.0 W05-W10 N21-N24
      7  363 1270 (N29E41 L=336)  1528    341        2 ( 264/1180- 453/1380) -38.7 E54-E33 N21-N37
      8 1170 1290 (N14W14 L= 31)   660    139        0 (1096/1255-1232/1332) -49.4 W10-W17 N11-N18
      9 1898 1121 (S10W67 L= 84)  1465    160        0 (1876/1040-1925/1196)  41.6 W62-W72 S14-S06
     10  532 1170 (N19E28 L=349)   245     47        0 ( 516/1135- 551/1192) -90.0 E29-E26 N17-N20
     11  198 1071 (N19E59 L=318)  2589    297        5 ( 119/ 991- 314/1179) -26.6 E73-E47 N13-N27
     12  420 1096 (N17E37 L=340)   783    126        1 ( 373/1058- 472/1134)   0.0 E41-E32 N15-N19
     13  451  805 (S00E40 L=337)   488     70        1 ( 411/ 781- 494/ 832)  -9.5 E42-E36 S02-N01
     14  529  761 (S04E36 L=342)   237     57        1 ( 512/ 729- 544/ 790)  36.9 E37-E33 S06-S02
     15 1055  675 (S19E05 L= 12)   330    149        1 (1037/ 612-1071/ 755)  58.0 E07-E02 S22-S14
     16 1303  634 (S27W10 L= 27)   295     45        1 (1281/ 601-1322/ 664) -90.0 W08-W11 S28-S24
     17  467  587 (S13E46 L=331)   436    148        0 ( 436/ 520- 497/ 649)  37.9 E50-E41 S17-S10
     18  514  528 (S18E45 L=332)   290     24        0 ( 494/ 512- 539/ 546) -45.0 E46-E43 S18-S16
     19  327  450 (S19E70 L=307)  1918     75        1 ( 286/ 412- 368/ 486) -16.7 E76-E64 S22-S16
  

    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