Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1716 1558 (N22W70 L=273)  1364    167        1 (1639/1531-1793/1580) -20.9 W61-W83 N17-N25
      2 1470 1528 (N27W43 L=246)   640    270        1 (1396/1484-1527/1566) -33.7 W37-W48 N24-N30
      3 1515 1443 (N21W43 L=245)   393     41        1 (1486/1421-1542/1471) -11.3 W39-W44 N19-N23
      4  879 1269 (N24E05 L=198)  1060    264        2 ( 798/1173- 940/1373) -70.3 E08-E01 N17-N31
      5 1382  824 (S12W18 L=221)   511    138        0 (1353/ 755-1426/ 875)  56.3 W15-W21 S15-S09
      6 1039  805 (S06E03 L=199)   440    178        0 (1002/ 745-1112/ 840)   9.5 E05-W01 S09-S04
      7  801  744 (S06E18 L=184)   595    192        2 ( 725/ 701- 890/ 788) -33.7 E22-E13 S10-S02
      8 1629  728 (S24W35 L=238)   364    128        0 (1578/ 689-1667/ 768)  14.0 W30-W38 S25-S21
      9 1769  711 (S29W50 L=252)   459     77        0 (1739/ 674-1802/ 748)  23.2 W46-W53 S30-S27
     10 1040  489 (S27E09 L=193)   534    207        0 ( 974/ 421-1118/ 557)  28.3 E15-E02 S30-S23
     11 1228  459 (S33W02 L=205)   439    122        0 (1188/ 418-1276/ 507)  35.5 E01-W06 S35-S30
     12  641  432 (S24E39 L=163)   282     55        1 ( 616/ 410- 662/ 462)  31.0 E41-E36 S25-S22
  

    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