Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1284 1685 (N40W40 L=171)   284     66        0 (1248/1672-1317/1700)  -8.1 W36-W43 N38-N41
      2 1638 1652 (N27W69 L=200)   517     62        1 (1612/1631-1672/1669) -26.6 W65-W73 N24-N28
      3  859 1557 (N43W02 L=133)   237     32        1 ( 836/1536- 882/1574) -33.7 W00-W03 N41-N44
      4  812 1477 (N39E04 L=127)   710    162        1 ( 791/1416- 853/1554)  63.4 E05-W00 N34-N43
      5 1046 1533 (N37W15 L=145)   267     76        1 (1014/1505-1086/1554) -11.3 W12-W17 N35-N37
      6  827 1200 (N21E08 L=123)  1198    341        1 ( 755/1073- 931/1335) -72.5 E11-E04 N11-N30
      7  667  927 (N08E23 L=108)   204     56        0 ( 632/ 912- 690/ 943)   0.0 E25-E21 N07-N09
      8  143  901 (N13E71 L= 60)   861    110        2 ( 115/ 859- 171/ 930) -68.2 E77-E66 N10-N15
      9  912  836 (S01E10 L=121)   279     93        1 ( 859/ 817- 954/ 854)   0.0 E13-E07 S02-N00
     10 1629  737 (S24W33 L=164)  1525    298        2 (1546/ 603-1695/ 847)  42.7 W25-W39 S30-S18
     11  490  706 (S02E40 L= 91)   381     65        1 ( 477/ 666- 510/ 748)  76.0 E41-E38 S04-N00
     12 1151  669 (S16E00 L=131)   681    109        6 (1114/ 623-1192/ 728)  82.9 E01-W02 S19-S11
     13  989  600 (S17E11 L=120)   231     72        0 ( 948/ 589-1021/ 612) -26.6 E13-E09 S18-S16
     14  872  537 (S19E20 L=111)   652    241        2 ( 817/ 482- 919/ 607) -83.7 E22-E16 S22-S13
     15 1662  549 (S39W37 L=167)   616    115        2 (1610/ 512-1703/ 589)  12.5 W31-W41 S40-S36
     16 1041  505 (S24E11 L=120)   285     96        0 (1009/ 468-1065/ 548)  45.0 E13-E08 S26-S21
     17  946  313 (S36E24 L=107)  5334    753       10 ( 668/ 167-1225/ 508)  13.0 E65-W00 S44-S28
     18  388  432 (S18E63 L= 68)   473     55        0 ( 364/ 411- 415/ 452)  11.3 E67-E57 S19-S16
     19 1579  357 (S54W32 L=163)   767     72        1 (1552/ 337-1612/ 385) -32.5 W27-W38 S57-S49
  

    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