Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1373 1776 (N53W49 L= 52)   905     89        0 (1332/1747-1412/1803)   7.1 W41-W58 N51-N55
      2 1265 1769 (N54W35 L= 38)   558     62        0 (1239/1739-1294/1792)  26.6 W30-W40 N51-N56
      3 1604 1671 (N41W66 L= 69)   656     86        1 (1571/1643-1646/1694) -39.8 W61-W70 N38-N43
      4  373 1535 (N41E58 L=305)   548     78        2 ( 327/1520- 412/1552)   4.8 E65-E52 N39-N41
      5  864 1451 (N33E09 L=354)  1154    243        2 ( 765/1391- 955/1532) -35.5 E16-E02 N28-N39
      6  680 1471 (N35E23 L=340)   315     87        0 ( 660/1424- 701/1514)  80.5 E24-E21 N32-N38
      7  920 1253 (N20E05 L=358)   210     88        0 ( 899/1220- 957/1274)  36.9 E06-E02 N17-N20
      8 1364 1052 (N04W22 L= 24)   263     76        0 (1336/1029-1406/1073)  21.8 W19-W24 N02-N05
      9 1744  851 (S12W50 L= 53)   690    108        0 (1710/ 808-1788/ 885)  33.7 W46-W54 S14-S10
     10  402  520 (S25E53 L=310)  1220    245        0 ( 332/ 448- 465/ 570) -45.0 E61-E48 S30-S21
     11  916  508 (S28E11 L=352)   409     95        0 ( 868/ 493- 974/ 524)   0.0 E14-E06 S28-S26
     12 1239  386 (S40W12 L= 15)   333     51        1 (1210/ 368-1270/ 407)  26.6 W09-W14 S41-S38
  

    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