Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  533 1763 (N25E62 L=237)   498     31        1 ( 515/1739- 553/1784) -45.0 E66-E59 N22-N26
      2 1087 1738 (N38E20 L=279)   595    121        0 (1040/1705-1140/1768)  -5.7 E25-E15 N36-N39
      3 1318 1623 (N38E00 L=299)   326     86        0 (1283/1598-1361/1643)   0.0 E02-W03 N35-N39
      4  576 1539 (N11E42 L=257)   925    204        1 ( 532/1468- 629/1641) -28.1 E50-E35 N07-N16
      5 1455 1334 (N23W16 L=316)   278     91        1 (1416/1302-1486/1367) -11.3 W13-W18 N21-N24
      6 1721 1220 (N25W38 L=338)   372    106        0 (1686/1176-1753/1242) -36.9 W35-W41 N23-N26
      7 1876 1097 (N24W57 L=356)   491    101        0 (1839/1054-1907/1138) -11.3 W51-W62 N22-N25
      8  668 1057 (S14E21 L=278)   371     71        1 ( 642/1024- 695/1093)  90.0 E22-E19 S15-S11
      9 1747  945 (N11W47 L=346)   659    128        1 (1722/ 885-1772/1005) -81.9 W44-W48 N07-N14
     10 1025  916 (S12W03 L=302)   476    185        1 ( 967/ 845-1083/ 987) -24.0 E01-W08 S14-S09
     11 1810  873 (N09W55 L=354)  2114    298        3 (1751/ 766-1873/ 984)  78.7 W50-W60 N02-N17
     12  934  829 (S19W00 L=300)   396    132        2 ( 871/ 802- 991/ 854)   0.0 E03-W04 S20-S17
     13  809  732 (S28E04 L=295)   469    152        1 ( 757/ 686- 872/ 761)  54.5 E07-W00 S32-S25
     14  210  690 (S47E63 L=236)   669     59        1 ( 176/ 670- 245/ 715)   7.1 E72-E56 S48-S45
     15 1418  620 (S17W36 L=335)   380     90        0 (1372/ 603-1462/ 639)  31.0 W32-W38 S18-S15
     16  529  547 (S49E21 L=278)   912    191        1 ( 449/ 506- 610/ 586)   0.0 E30-E13 S50-S46
     17  851  462 (S44W08 L=307)   355     90        1 ( 812/ 436- 893/ 493)   8.1 W04-W11 S45-S41
     18 1577  332 (S26W69 L=  8)   489     52        0 (1547/ 310-1601/ 351)   5.2 W63-W74 S26-S24
     19 1025  166 (S58W49 L=348)   465     48        0 (1007/ 142-1046/ 199) -17.1 W42-W55 S59-S55
  

    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