Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  826 1830 (N44E46 L=344)  1824    208        4 ( 733/1776- 903/1873)  -2.5 E58-E34 N41-N48
      2 1556 1618 (N51W28 L= 58)  3248    528        6 (1329/1492-1707/1750) -10.3 W01-W47 N44-N55
      3 1034 1577 (N31E14 L= 17)   617    146        0 ( 982/1519-1079/1619) -24.0 E18-E09 N27-N33
      4  782 1521 (N21E28 L=  3)   586    214        0 ( 758/1418- 820/1603) -72.9 E30-E24 N14-N27
      5  638 1449 (N13E34 L=356)   698    134        2 ( 619/1390- 669/1521) -45.0 E38-E30 N10-N17
      6  452 1412 (N07E47 L=344)   339     62        0 ( 430/1383- 482/1435)  76.0 E48-E44 N05-N09
      7 1056 1309 (N14E05 L= 26)   317    102        1 (1022/1291-1119/1345)  18.4 E06-W00 N12-N16
      8 1785 1242 (N29W50 L= 81)   883    142        1 (1738/1201-1842/1290) -13.0 W44-W57 N27-N31
      9  173 1143 (S14E67 L=324)   884     93        2 ( 134/1114- 222/1163)   8.7 E73-E60 S15-S12
     10  263  969 (S22E53 L=338)  2214    425        4 ( 223/ 833- 338/1148) -81.0 E57-E44 S30-S10
     11 1880 1075 (N22W63 L= 94)   714    132        0 (1842/1039-1933/1124)  -6.7 W57-W74 N20-N23
     12  879  926 (S12E06 L= 25)   298    106        0 ( 831/ 909- 940/ 940)  26.6 E08-E02 S13-S09
     13  465  786 (S30E32 L=359)   341    118        0 ( 424/ 744- 512/ 816)  40.6 E35-E28 S32-S26
     14  769  703 (S27E08 L= 22)   267    165        0 ( 695/ 688- 853/ 713)  24.4 E13-E02 S29-S24
     15 1183  468 (S32W26 L= 57)  1627    326        2 (1024/ 435-1330/ 521)  29.1 W15-W34 S37-S26
     16  667  400 (S52E10 L= 20)   327     48        0 ( 635/ 388- 693/ 417)  38.7 E13-E07 S54-S50
     17  690  252 (S67E04 L= 27)  1877    236        4 ( 590/ 219- 804/ 287)   8.4 E22-W12 S70-S61
  

    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