Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1201 1893 (N56W51 L=348)   705     64        0 (1179/1856-1217/1923)  25.0 W43-W60 N52-N59
      2 1212 1693 (N39W29 L=327)  1405    233        2 (1092/1658-1315/1727) -13.2 W19-W36 N35-N41
      3 1583 1574 (N24W52 L=350)   363     68        0 (1547/1547-1611/1596) -45.0 W49-W53 N21-N25
      4 1475 1531 (N23W41 L=338)   664     77        1 (1448/1495-1514/1568)  26.6 W37-W44 N21-N25
      5  907 1460 (N29W01 L=298)  1538    512        2 ( 666/1432-1142/1482) -17.9 E16-W15 N24-N34
      6 1355 1391 (N17W27 L=325)   567    174        0 (1270/1371-1437/1409) -15.3 W21-W32 N14-N18
      7  587 1367 (N29E23 L=275)   221     71        0 ( 553/1340- 612/1388) -36.9 E25-E21 N27-N30
      8  352 1331 (N31E43 L=254)  1208    239        2 ( 246/1286- 457/1381) -16.9 E56-E33 N27-N35
      9  943 1349 (N22W01 L=298)   515    213        1 ( 835/1328-1031/1366) -26.6 E05-W05 N19-N24
     10  920 1240 (N15E03 L=295)   892    384        1 ( 789/1133-1062/1324) -49.1 E09-W04 N06-N22
     11 1845 1069 (S10W58 L=356)   626    128        0 (1832/ 999-1858/1137)  69.4 W56-W60 S14-S06
     12  731  958 (N02E18 L=279)   295    111        1 ( 679/ 921- 771/ 983)  15.9 E22-E15 N00-N03
     13 1574  921 (S14W33 L=330)   270    121        0 (1524/ 888-1623/ 960)  20.6 W28-W36 S15-S12
     14 1190  753 (S18W05 L=303)   235     76        0 (1164/ 725-1230/ 784)  11.3 W03-W08 S19-S16
     15  932  731 (S14E11 L=287)   361    121        0 ( 885/ 693- 970/ 773) -63.4 E12-E09 S17-S11
     16 1229  744 (S19W08 L=305)   206     47        0 (1207/ 722-1257/ 763)  14.0 W05-W09 S20-S17
     17 1023  493 (S31E10 L=287)  1091    233        1 ( 990/ 386-1069/ 609)  54.5 E15-E04 S38-S24
  

    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