Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  776 1818 (N57E40 L=308)  1521    119        3 ( 736/1756- 826/1854) -68.2 E48-E34 N50-N61
      2  665 1821 (N55E57 L=292)   587     52        0 ( 636/1806- 698/1835)   0.0 E63-E50 N54-N56
      3  955 1596 (N39E10 L=338)   295     67        1 ( 932/1573- 982/1612)  45.0 E12-E07 N37-N40
      4  507 1494 (N28E43 L=305)   670    124        1 ( 438/1461- 559/1543) -12.1 E52-E38 N25-N30
      5 1492 1413 (N30W32 L= 20)   632    210        0 (1441/1340-1561/1482) -45.0 W28-W36 N25-N34
      6  313 1401 (N20E59 L=290)  1593    171        4 ( 260/1329- 359/1478) -27.9 E69-E52 N15-N24
      7  774 1322 (N19E18 L=330)   430    168        0 ( 709/1265- 829/1360) -29.1 E23-E14 N15-N20
      8 1338 1078 (N07W19 L=  8)   437     84        1 (1306/1039-1367/1114)  63.4 W17-W21 N04-N09
      9 1249  698 (S17W17 L=  5)   282     68        1 (1225/ 658-1273/ 733)  78.7 W15-W18 S19-S14
     10 1444  357 (S40W43 L= 32)   459     46        1 (1427/ 334-1467/ 385) -29.7 W39-W46 S42-S37
     11 1091  319 (S46W12 L=  1)   331     48        0 (1066/ 301-1116/ 340) -31.0 W09-W14 S47-S44
     12 1253  298 (S47W28 L= 17)   404     21        0 (1233/ 273-1274/ 321) -33.7 W25-W31 S48-S44
  

    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