Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1431 1785 (N38W64 L=303)   960    138        1 (1383/1749-1485/1807)  -3.2 W55-W73 N35-N39
      2 1706 1484 (N13W61 L=300)  1504    224        2 (1592/1444-1802/1540) -34.5 W52-W70 N07-N19
      3 1024 1347 (N25W09 L=248)   547    150        1 ( 951/1302-1086/1398) -58.0 W05-W11 N21-N29
      4 1885 1280 (S04W70 L=309)   562     79        1 (1866/1241-1898/1323)  18.4 W65-W74 S06-S02
      5 1052 1255 (N19W07 L=246)   437    175        2 ( 970/1219-1131/1277) -11.3 W02-W12 N17-N20
      6 1606 1144 (S02W38 L=276)   713    137        0 (1583/1061-1632/1205)  50.2 W34-W40 S05-N01
      7 1068 1071 (N09W04 L=242)   259    131        0 (1034/1016-1101/1134)  39.8 W00-W07 N06-N11
      8  321 1064 (N24E47 L=192)   327     37        2 ( 294/1039- 344/1091) -76.0 E49-E44 N22-N26
      9  260  905 (N16E55 L=184)   581    108        0 ( 225/ 856- 289/ 941) -63.4 E58-E52 N12-N18
     10  190  855 (N14E64 L=174)  1050    103        1 ( 169/ 804- 212/ 910)  50.2 E67-E61 N11-N17
     11 1856  869 (S26W55 L=294)   758    111        0 (1811/ 849-1895/ 901) -18.4 W49-W61 S28-S23
     12  463  838 (N09E39 L=200)   400    101        0 ( 434/ 793- 495/ 879) -71.6 E40-E37 N05-N11
     13  236  788 (N10E61 L=178)   475     51        0 ( 208/ 772- 274/ 804) -15.9 E63-E56 N08-N10
     14  433  535 (S08E52 L=187)   781    161        0 ( 381/ 470- 460/ 598) -83.7 E55-E49 S12-S03
     15 1394  398 (S42W07 L=246)  3218    670        6 (1134/ 241-1644/ 570)   9.3 E18-W31 S47-S37
     16 1182  516 (S27E04 L=235)   393     85        0 (1162/ 479-1207/ 559) -90.0 E05-E02 S29-S23
     17  630  476 (S16E39 L=200)  1306    234        6 ( 570/ 399- 694/ 553)  73.3 E44-E34 S21-S11
     18  777  427 (S22E32 L=207)   533    102        0 ( 733/ 399- 805/ 469)  15.9 E35-E28 S23-S19
  

    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