Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  704 1437 (N08E30 L=284)   980    303        0 ( 670/1295- 733/1574) -60.6 E35-E25 N00-N16
      2 1808 1431 (N45W47 L=  1)   541    117        0 (1791/1399-1821/1466)  54.5 W44-W49 N41-N48
      3 1678 1338 (N32W33 L=347)   378     75        0 (1646/1302-1706/1366) -18.4 W29-W35 N30-N33
      4  881 1298 (N05E15 L=299)   240     77        0 ( 838/1274- 921/1315)  36.9 E17-E12 N02-N06
      5  401 1152 (S15E42 L=273)  1766    437        2 ( 235/1103- 543/1177)  16.3 E56-E31 S19-S11
      6 1874  748 (N06W69 L= 24)   801    133        1 (1845/ 731-1907/ 769)   9.5 W64-W76 N04-N07
      7 1515  692 (S10W39 L=353)   383    109        0 (1493/ 637-1542/ 746) -81.9 W37-W40 S13-S06
      8 1217  503 (S29W29 L=343)   738    132        2 (1145/ 466-1264/ 552)  69.4 W26-W31 S33-S25
      9 1686  508 (S13W62 L= 16)   618     70        1 (1655/ 479-1723/ 542)  59.0 W58-W64 S15-S10
     10 1085  465 (S36W23 L=337)   239     28        0 (1068/ 444-1103/ 491) -45.0 W21-W24 S37-S34
     11 1520  350 (S27W58 L= 12)  2037    381        3 (1438/ 296-1595/ 395)   3.0 W49-W68 S28-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