Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  483 1560 (N53E37 L= 73)   557     63        1 ( 464/1519- 505/1590) -56.3 E40-E33 N49-N55
      2  718 1561 (N47E11 L= 99)   698    113        1 ( 663/1535- 769/1581) -29.1 E15-E06 N44-N49
      3 1788 1451 (N12W70 L=180)   791     81        0 (1770/1408-1803/1486)  21.8 W67-W77 N09-N13
      4  671 1275 (N29E18 L= 92)   237     77        0 ( 651/1240- 700/1302)  36.9 E19-E15 N26-N30
      5  314 1130 (N25E49 L= 61)   872    217        1 ( 243/1051- 366/1203) -42.5 E57-E44 N19-N30
      6  180 1114 (N25E66 L= 44)  1721    316        1 ( 118/1000- 239/1197)  33.7 E77-E58 N18-N30
      7  790  818 (S01E18 L= 92)   651    196        1 ( 742/ 728- 843/ 894) -90.0 E19-E16 S06-N04
      8 1056  736 (S11E04 L=106)   256     88        1 (1024/ 696-1086/ 767)  45.0 E06-E01 S12-S08
      9 1709  712 (S28W42 L=152)   438     41        1 (1687/ 685-1729/ 743)  90.0 W39-W43 S30-S26
     10  404  417 (S21E64 L= 46)   841     75        1 ( 371/ 384- 433/ 440) -76.0 E67-E59 S23-S18
     11  489  368 (S25E58 L= 52)  1017    106        0 ( 460/ 319- 528/ 415)  69.4 E61-E54 S29-S21
  

    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