Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  446 1602 (N19E58 L=111)  3399    399        8 ( 342/1520- 544/1684)  52.4 E67-E50 N13-N26
      2 1647 1647 (N58W45 L=213)   683     65        1 (1631/1610-1671/1673)  31.6 W39-W55 N53-N61
      3  363 1504 (N11E59 L=109)   987    141        1 ( 298/1466- 411/1546)  29.1 E66-E55 N08-N14
      4 1582 1495 (N41W29 L=198)   342     54        1 (1563/1472-1606/1521) -33.7 W27-W31 N39-N42
      5  581 1488 (N15E40 L=128)   282     72        1 ( 561/1453- 611/1511)  90.0 E41-E38 N12-N16
      6 1695 1340 (N33W40 L=209)   481     70        1 (1671/1309-1718/1377) -45.0 W37-W42 N31-N35
      7  175 1238 (S08E70 L= 99)   684     68        0 ( 140/1212- 209/1261)  -4.4 E77-E63 S08-S06
      8  150  929 (S27E69 L=100)   735     77        1 ( 130/ 894- 176/ 969)  39.8 E74-E64 S29-S24
      9  868  721 (S24E02 L=166)   282     94        1 ( 846/ 679- 895/ 760) -78.7 E03-W00 S26-S21
     10  781  669 (S29E07 L=162)   222     77        0 ( 739/ 648- 816/ 690)  -9.5 E09-E03 S29-S28
     11 1112  631 (S23W15 L=184)   246     76        1 (1071/ 610-1144/ 652)  11.3 W12-W17 S24-S22
     12 1439  626 (S16W36 L=205)   248     58        0 (1410/ 600-1466/ 648)  63.4 W34-W37 S17-S13
     13 1309  618 (S19W28 L=196)   246     51        0 (1293/ 587-1321/ 646) -45.0 W26-W29 S21-S17
     14  322  586 (S47E51 L=118)  1251    145        2 ( 266/ 561- 390/ 628)   0.0 E59-E40 S49-S44
     15  760  538 (S38E05 L=163)   291     48        0 ( 737/ 510- 784/ 562)  76.0 E06-E03 S40-S36
     16 1553  372 (S27W62 L=230)  1346    163        2 (1502/ 328-1616/ 428)  83.7 W56-W65 S31-S22
     17 1157  392 (S37W29 L=197)   455     82        1 (1114/ 369-1191/ 412)   0.0 W24-W32 S38-S35
     18  779  376 (S50W01 L=169)   681    136        1 ( 754/ 317- 817/ 410) -40.6 E02-W05 S54-S46
     19  794  278 (S59W07 L=176)   927    111        1 ( 747/ 245- 851/ 314)  29.1 W01-W10 S62-S53
  

    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