Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1006 1839 (N62W34 L=281)   403     67        1 ( 978/1816-1041/1856)   0.0 W29-W40 N59-N63
      2 1451 1703 (N37W56 L=303)   563     67        1 (1416/1681-1485/1719) -15.9 W52-W60 N34-N38
      3 1520 1665 (N32W59 L=306)   473     39        1 (1496/1644-1551/1683)  -9.5 W55-W62 N31-N33
      4  555 1639 (N58E30 L=217)   675    101        1 ( 517/1597- 595/1677)  36.9 E34-E24 N54-N60
      5 1286 1479 (N28W28 L=275)   333     73        1 (1258/1453-1323/1509)   0.0 W25-W31 N26-N29
      6  351 1367 (N40E47 L=199)   779    114        1 ( 310/1339- 407/1391)   0.0 E52-E40 N38-N42
      7 1760  834 (S22W47 L=294)   502     79        2 (1736/ 798-1785/ 875)  53.1 W44-W49 S23-S19
      8 1717  644 (S34W43 L=290)   437     63        1 (1682/ 615-1746/ 676)   0.0 W39-W46 S35-S31
      9  479  565 (S11E46 L=200)   494     76        0 ( 437/ 540- 517/ 591)   7.1 E50-E42 S12-S09
     10  612  546 (S14E37 L=210)   400    114        1 ( 590/ 500- 638/ 587)  80.5 E38-E34 S17-S11
     11  579  449 (S20E44 L=203)   839    211        3 ( 512/ 415- 664/ 484) -14.0 E48-E36 S23-S17
  

    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