Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1054 1566 (N33W11 L=293)   240    132        0 (1025/1521-1090/1629) -83.7 W10-W12 N29-N38
      2 1893  969 (S16W64 L=346)   525     41        1 (1872/ 946-1911/ 996)  33.7 W61-W67 S16-S13
      3 1780  935 (S16W51 L=333)   315     35        1 (1762/ 910-1797/ 961) -71.6 W48-W52 S17-S14
      4  698  543 (S24E29 L=253)   256     46        1 ( 673/ 518- 719/ 570)  76.0 E30-E27 S25-S21
      5 1673  503 (S43W47 L=329)   592     61        1 (1650/ 472-1691/ 542) -59.0 W43-W49 S46-S40
      6  695  469 (S28E32 L=250)   380     97        1 ( 663/ 418- 722/ 505)  80.5 E34-E29 S31-S25
      7  993  466 (S33E12 L=270)   518     70        2 ( 967/ 420-1022/ 500)  45.0 E14-E09 S36-S30
      8 1303  445 (S41W11 L=293)   250     64        0 (1281/ 420-1317/ 480)  63.4 W08-W11 S42-S37
      9 1157  320 (S48E03 L=279)   559    101        0 (1117/ 302-1203/ 340)   0.0 E07-W00 S49-S45
  

    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