Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  493 1235 (N31E28 L= 37)  1870    343        1 ( 353/1182- 646/1275) -14.0 E41-E17 N25-N34
      2  420 1017 (N20E38 L= 28)  1234    306        3 ( 303/ 958- 551/1094) -37.4 E47-E29 N13-N26
      3  706  936 (N09E20 L= 45)  1363    375        3 ( 643/ 864- 780/1038)  84.3 E23-E14 N04-N14
      4  197  931 (N19E61 L=  4)  1363    225        2 ( 147/ 834- 274/1016) -47.3 E67-E54 N12-N25
      5 1823  896 (S26W51 L=116)   408     48        1 (1798/ 871-1842/ 925)  14.0 W47-W53 S26-S24
      6 1479  817 (S19W21 L= 86)   302    128        0 (1459/ 752-1492/ 874)  56.3 W18-W22 S21-S15
      7  856  772 (S04E16 L= 49)   222     64        0 ( 821/ 750- 877/ 798)  14.0 E18-E14 S05-S02
      8 1750  600 (S43W43 L=108)   625     72        2 (1719/ 570-1785/ 629)   7.1 W39-W47 S44-S40
      9 1613  526 (S43W27 L= 92)   361     75        0 (1583/ 508-1647/ 546)   0.0 W23-W30 S43-S41
     10  451  373 (S19E62 L=  3)  2913    350        6 ( 327/ 255- 544/ 497)  87.1 E71-E56 S28-S08
     11  922  281 (S37E33 L= 33)  1830    231        3 ( 834/ 200- 999/ 360)  20.0 E44-E22 S41-S32
     12  777  197 (S39E52 L= 14)   419     75        0 ( 763/ 176- 804/ 233)  50.2 E55-E47 S41-S35
  

    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