Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  378 1601 (N30E66 L= 28)   570     59        1 ( 360/1571- 394/1636) -33.7 E69-E62 N28-N32
      2 1408 1599 (N44W25 L=118)   453     64        1 (1368/1581-1444/1618)  20.6 W20-W28 N42-N45
      3 1203 1364 (N24W08 L=101)  1156    486        1 (1067/1226-1342/1497) -35.0 E02-W18 N16-N31
      4  316 1442 (N19E61 L= 33)   710    144        1 ( 272/1394- 377/1474) -13.2 E68-E51 N16-N20
      5 1648 1414 (N32W44 L=138)   534    133        0 (1598/1372-1688/1460) -35.5 W40-W47 N29-N35
      6 1309 1302 (N21W16 L=109)   439    103        2 (1270/1280-1361/1337) -15.9 W12-W19 N19-N23
      7 1128 1239 (N15W04 L= 98)   286    128        0 (1113/1182-1142/1298) -76.0 W02-W05 N11-N19
      8  641 1000 (S04E24 L= 70)   246     75        0 ( 604/ 982- 669/1026)  26.6 E26-E21 S05-S02
      9 1046  720 (S17W05 L= 98)   550    141        3 ( 991/ 683-1113/ 759)  15.9 W01-W08 S20-S14
     10 1515  607 (S20W39 L=132)   797    154        4 (1448/ 551-1595/ 632)  14.0 W33-W45 S23-S17
     11  336  553 (S38E55 L= 39)  1201    152        1 ( 278/ 511- 400/ 594)  15.3 E67-E45 S41-S33
     12 1546  369 (S36W56 L=149)   431     68        1 (1520/ 351-1582/ 392)  31.0 W52-W58 S37-S34
     13 1051  196 (S60W20 L=114)   833    121        0 (1027/ 151-1068/ 257) -60.3 W15-W26 S66-S52
  

    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