Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  370 1571 (N28E64 L= 30)  1073    146        1 ( 343/1503- 396/1637) -48.4 E68-E59 N23-N32
      2 1410 1602 (N44W25 L=118)   476     58        1 (1379/1582-1445/1620)  23.2 W21-W28 N42-N45
      3 1212 1360 (N24W09 L=102)  1116    453        0 (1083/1226-1342/1481) -34.4 E01-W18 N16-N31
      4  297 1452 (N19E64 L= 29)   532     80        1 ( 275/1422- 331/1472)  21.8 E68-E58 N17-N20
      5 1652 1414 (N32W45 L=138)   543    121        0 (1601/1372-1689/1460) -40.6 W40-W48 N29-N35
      6 1309 1301 (N21W16 L=109)   503    101        2 (1273/1280-1370/1336)  -9.5 W13-W19 N19-N23
      7 1133 1238 (N15W05 L= 98)   256    110        0 (1117/1184-1146/1293) -80.5 W03-W05 N12-N18
      8 1053  720 (S17W05 L= 98)   529    146        1 (1001/ 683-1119/ 758)  33.7 W01-W09 S19-S14
      9 1516  612 (S19W39 L=132)   641    181        2 (1450/ 588-1596/ 630)  10.3 W33-W45 S20-S17
     10  338  552 (S38E54 L= 39)  1236    154        1 ( 279/ 511- 403/ 594)  16.3 E68-E44 S41-S33
     11 1054  195 (S60W21 L=114)   827    120        0 (1028/ 148-1071/ 256) -57.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