Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1097 1859 (N51E24 L= 82)   760     84        1 (1046/1841-1144/1875)  39.8 E27-E19 N47-N53
      2  671 1780 (N32E49 L= 57)   417     43        0 ( 649/1749- 691/1806) -18.4 E53-E44 N29-N34
      3 1441 1690 (N47W12 L=118)   517     54        1 (1412/1654-1474/1713)   9.5 W09-W15 N43-N49
      4  645 1651 (N22E42 L= 64)   522    104        0 ( 602/1615- 691/1676)  45.0 E43-E38 N19-N24
      5  956 1560 (N23E17 L= 88)   312    142        0 ( 924/1489- 989/1607) -40.6 E20-E13 N18-N25
      6  510 1476 (N09E43 L= 63)   295     86        1 ( 483/1446- 545/1514)  68.2 E44-E41 N06-N11
      7  142 1117 (S17E70 L= 36)   586     90        1 ( 123/1078- 162/1162)  51.3 E72-E66 S19-S14
      8  135  970 (S26E70 L= 35)   567     73        3 ( 116/ 940- 159/1004)  63.4 E74-E66 S28-S24
      9  292  606 (S48E56 L= 50)   373     39        0 ( 273/ 589- 321/ 626)   7.1 E59-E50 S49-S46
     10 1480  434 (S28W50 L=155)   963    230        1 (1420/ 358-1577/ 498)  79.7 W47-W53 S33-S22
     11  904  267 (S58W20 L=126)   365     62        1 ( 880/ 243- 931/ 297) -14.0 W16-W24 S60-S56
     12 1123  158 (S58W55 L=161)   987    109        0 (1095/ 120-1159/ 197) -25.0 W46-W66 S61-S54
  

    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