Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1684 1432 (N26W55 L=346)   932    161        0 (1632/1366-1723/1498) -77.5 W52-W56 N21-N30
      2 1107 1423 (N29W07 L=299)   616    140        2 (1082/1362-1133/1487)  76.0 W05-W10 N25-N33
      3  283 1416 (N30E63 L=228)   643    131        0 ( 247/1371- 337/1456)  36.9 E67-E57 N26-N32
      4  292 1340 (N25E58 L=234)   532     63        0 ( 273/1305- 313/1373)  78.7 E59-E56 N22-N27
      5 1184 1328 (N23W12 L=303)   494     99        0 (1133/1309-1237/1348) -15.9 W08-W15 N21-N24
      6  949  867 (S05E05 L=286)   202     42        1 ( 929/ 841- 974/ 885)  33.7 E06-E03 S06-S03
      7 1747  853 (S11W51 L=343)   544     17        1 (1722/ 825-1772/ 877)  71.6 W48-W53 S12-S08
      8 1062  806 (S09W02 L=293)   390     86        2 (1032/ 771-1091/ 840)  53.1 W00-W03 S11-S07
      9  874  783 (S10E10 L=281)   330     93        0 ( 853/ 745- 897/ 830) -90.0 E11-E08 S12-S07
     10 1331  744 (S14W19 L=310)   257     74        0 (1314/ 709-1342/ 784)  90.0 W17-W19 S16-S11
     11 1694  606 (S27W51 L=342)   879    120        1 (1664/ 579-1738/ 641)   7.1 W47-W55 S28-S24
     12 1689  496 (S35W56 L=348)   447     58        0 (1655/ 479-1712/ 512)   0.0 W51-W60 S36-S34
  

    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