Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1259 1605 (N40W01 L=218)  1131    209        2 (1189/1524-1300/1679) -32.5 E06-W05 N34-N44
      2  706 1472 (N17E31 L=186)  1334    510        4 ( 552/1345- 865/1634) -16.9 E50-E16 N12-N24
      3  740 1559 (N23E32 L=185)   326    107        0 ( 709/1523- 775/1612) -26.6 E36-E28 N21-N26
      4 1176 1451 (N27E01 L=216)  1000    406        1 (1071/1287-1233/1596) -43.3 E11-W06 N17-N34
      5 1085 1372 (N19E04 L=213)   346    145        1 (1063/1305-1103/1438) -71.6 E06-E02 N14-N23
      6 1534 1289 (N25W27 L=244)   640    176        1 (1507/1230-1592/1331)  54.5 W24-W31 N21-N28
      7  234 1243 (S06E62 L=156)   443     45        1 ( 210/1225- 254/1262)  -8.1 E65-E58 S07-S04
      8  222 1111 (S14E59 L=158)  2115    305        4 ( 180/ 974- 264/1215) -57.3 E67-E54 S22-S08
      9 1782 1063 (N18W51 L=268)   860    147        2 (1762/ 984-1801/1137) -84.3 W48-W53 N12-N22
     10  758  964 (S12E14 L=203)   407    145        1 ( 712/ 932- 801/1006)  -9.5 E17-E11 S13-S09
     11 1187  911 (S05W12 L=229)   388    168        1 (1117/ 854-1220/ 959) -26.6 W07-W15 S07-S03
     12  900  761 (S21E01 L=216)   552    262        1 ( 788/ 723-1001/ 784)  39.3 E07-W04 S25-S16
     13  331  701 (S38E46 L=171)  3236    387        5 ( 166/ 663- 483/ 754)   8.3 E85-E30 S43-S32
     14 1615  396 (S25W66 L=283)   531     61        1 (1602/ 370-1631/ 432) -26.6 W61-W71 S26-S22
     15  695  399 (S51E08 L=209)   564    114        1 ( 647/ 378- 751/ 429)   0.0 E13-E02 S52-S49
     16 1570  365 (S28W65 L=282)   539     49        0 (1548/ 348-1600/ 386)   0.0 W60-W69 S28-S26
     17 1498  335 (S32W60 L=277)   737     51        3 (1471/ 308-1527/ 366) -20.6 W55-W64 S33-S30
  

    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